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I. OVERVIEW 

These lectures are about some recent developments concerning the physics of the strong 
interaction, Quantumchromodynamics (QCD), in extremes of temperature and baryon den- 
sity. Here "extremes" means temperatures of the order of 10^^ K or 100 MeV and densities 
of about a few times nuclear matter density no = 0.153 fm^ = (105 MeV)^. [| 

The following schematic phase diagram gives an idea about what the behavior of QCD in 
thermal and chemical equilibrium may look like as a function of temperature and chemical 
potential of quark number density. We will use it here to draw the attention to some aspects 
which will be explained in more detail in these lectures and to point out their experimental 
relevance. 
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Shortly after the discovery of asymptotic freedom |l[], it was realized [|] that at sufficiently 
high temperature or density QCD may differ in important aspects from the corresponding 
zero temperature or vacuum properties. If one considers the phase structure of QCD as 
a function of temperature one expects around a critical temperature of Tc ~ 150 MeV two 
qualitative changes. First, the quarks and gluons which are confined into hadrons at zero 
temperature can be excited independently at sufficiently high temperature. The hot state 
is conventionally called the quark— gluon plasma. The second aspect has to do with the 
fact that the vacuum of QCD contains a condensate of quark-ant iquark pairs, (ipip) ^ 0, 
which spontaneously breaks the (approximate) chiral symmetry of QCD and has profound 



^We will work in units where 7i = c = 1 such that the mass of a particle is equal to its rest energy 
{mc^) and also to its inverse Compton wavelength (mc/h). It is sometimes convenient to express 
length in the unit of a fermi (fm = 10^^"^ cm) because it is the order of the dimension of a nucleon. 
The conversion to units of energy is easily done through (1 MeV)~^hc = 197.33 fm . 



implications for the hadron spectrum. At high temperature this condensate is expected to 
melt, i.e. (^^) — 0, which signals the chiral phase transition. 

QCD at nonzero baryon density is expected to have a rich phase structure with different 
possible phase transitions as the density varies. Since nucleons as bound states of quarks 
have a characteristic size tn — Ifm ~ (200MeV)~^, for very high baryon density there 
is not enough space to form nucleons and one expects a new quark matter phase. In 
this phase, similar to the case of very high temperature, confinement is not expected to 
play an important role and the quark-antiquark condensate is absent. However, the high 
density state is far from trivial. It has been realized early ^ that at very high density, 
where perturbative QCD can be applied, quark matter behaves as a color superconductor: 
Cooper pairs of quarks condense, opening up an energy gap at the quark Fermi surface. 
Recent investigations 0] at intermediate densities using effective models show the formation 
of quark-quark condensates with phenomenologically significant gaps of order 100 MeV. 
A first study which takes into account the formation of condensates in the conventional 
quark-antiquark channel and in a superconducting quark-quark channel reveals a phase 
diagram similar to the one sketched in the above figure. New symmetry breaking schemes 
like color-flavor locking may also be relevant for the study of nuclear matter, if the 
latter is continuously connected to the quark matter phase [0. The nuclear matter state 
is intermediate between a gas of nucleons and quark matter and can be associated with a 
liquid of nucleons. 

Where do we encounter QCD at high temperatures and/or densities? The astrophysics 
of neutron stars provides a good testing ground for the exploration of very dense matter. 
Neutron stars are cold on the nuclear scale with temperatures of about 10^ to 10^ K (1 MeV = 
1.1065 X 10^°K). The range of densities is enormous. At the edge, where the pressure is 
zero, the density is that of ordinary iron. In contrast, at the center it may be a few times 
nuclear matter density. 

According to the standard hot big bang cosmology the high temperature transition must 
have occured during the evolution of the early universe. For most of its evolution the 
early universe was to a good approximation in thermal equilibrium and the transition took 
place about a microsecond after the big bang, where the temperature dropped to the order 
of 100 MeV. 

A promising prospect is to reproduce QCD phase transitions in heavy— ion collisions 
in the laboratory. Large efforts in ultra-relativistic heavy ion collision experiments focus 
on the observation of signatures for a high temperature and/or density transition |Q]. In- 
tensive searches have been performed at the AGS accelerator (BNL) and at the CERN 
SPS accelerator and soon also at the RHIC collider and the future LHC. In a relativistic 
heavy-ion collision, one may create a region of the high temperature or density phase in 
approximate local equilibrium. Depending on the initial density and temperature, when 
this region expands and cools it will traverse the phase transition at different points in the 
phase diagram. A main challenge is to find distinctive, qualitative signatures that such a 
transition has indeed occured during the course of the collision. 

Certain signatures rely on the observation that near a phase transition long-range cor- 
relations can occur. We are familiar with this phenomenon from condensed matter systems. 
An example is a ferromagnet which is heated above a certain temperature Tc where its 
magnetization disappears. The phase transition is second order and in the vicinity of T^ 



the correlation length between spins grows very large. Points in the phase diagram which 
correspond to second order phase transitions are called critical points. It can be argued 
p|^Tl| that in the theoretical limit of two massless quark flavors there is a line of critical 
points in the phase diagram of QCD (cf. the dashed line in the above phase diagram). In- 
deed, in the real world there are two quarks, the u and the d, which are particularly light. 
A large correlation length may then be responsible ||TO| , p!T| for the creation of large domains, 



in which the pion fleld has a non-zero expectation value pointing in a flxed direction in 
isospin space ("disoriented chiral condensate"). The observational consequences would be 
strong fluctuations in the number ratio of neutral to charged pions |TD|JTT[] . The light quarks, 
however, are not massless and it is a quantitative question if the correlation length is long 
enough to allow for distinctive signatures in a heavy-ion collision. 

Most strikingly, there is the possibility of a truely inflnite correlation length for a particu- 
lar temperature and density even for realistic quark masses. The corresponding critical point 
in the phase diagram marks the endpoint of a line of flrst order phase transitions. A familiar 
analogy in condensed matter physics comes from the boiling of water. For the liquid-gas 
system there exists an endpoint of flrst order transitions in the phase diagram where the 
liquid and the gaseous phase become indistinguishable and which exhibits the well-known 
phenomenon of critical opalescence. It is the precise analogue to critical opalescence in QCD 
at nonzero temperature and density which has drawn much attention recently [^,12|,13|. In 



QCD the corresponding critical point marks the endpoint of a line of flrst order transitions 
in which the quark-antiquark condensate (ip'^') drops discontinuously as the temperature or 
density is increased. In the theoretical limit of vanishing quark masses this critical endpoint 
becomes a tricritical point (cf. the above flgure). We note that a critical endpoint is also 
known for the nuclear gas-liquid transition for a temperature of about 10 MeV. Signatures 
and critical properties of this point have been studied through measurements of the yields 



of nuclear fragments in low energy heavy ion collisions ||14| , p!5 |. 

Apart from the phenomenological implications, a large correlation length near a critical 
point opens the possibility that the QCD phase transition is characterized by universal prop- 
erties. The notion of universality for critical phenomena is well-established in statistical 
physics. Universal properties are independent of the details, like short distance couphngs, 
of the model under investigation. They only depend on the symmetries, the dimensionality 
of space and the fleld content. As a consequence a whole class of models is described by the 
same universal scaling form of the equation of state in the vicinity of the critical point. The 
range of applicability typically covers very different physical systems in condensed matter 
physics and high temperature quantum fleld theory. The analogy between critical opales- 
cence in a liquid-gas system and QCD near a critical endpoint is indeed quantitatively 
correct: Both systems are in the universality class of the well-known Ising model. 

The thermodynamics described above is difficult to tackle analytically. A major problem 
is that for the relevant length scales the running QCD gauge coupling a^ is large making 
a perturbative approach unreliable. The universal QCD properties near critical points may 
be computed within a much simpler model in the same universality class. However, effective 
couplings near critical points are typically not small and the physics is characterized by 
nonanalytic behavior. Important information near the phase transition is also nonuniversal, 
like the critical temperature T^ or the overall size of a correlation length. The question how 



small rriu and ma would have to be in order to see a large correlation length near Tc at 
low density, and if it is realized for realistic quark masses, will require both universal and 
nonuniversal information. 

A very promising approach to treat these questions is the use of the Wilsonian 'exact 
renormalization group' applied to an effective field theory. More precisely, we employ 
an exact nonperturbative flow equation for a scale dependent effective action r^ ||16|| , 
which is the generating functional of the IP/ Green functions in the presence of an infrared 
momentum cutoff ~ k. In the language of statistical physics, F^ is a coarse grained free 
energy with a coarse graining length scale ~ k~^. The renormalization group flow for the 
average action F^ interpolates between a given short distance or classical action S and the 
standard effective action F, which is obtained by following the flow for Tk to k = 0. We 
will investigate in these lectures the QCD phase diagram within the linear quark meson 
model for two quark flavors. Truncated nonperturbative flow equations are derived at 
nonzero temperature and chemical potential. Whereas the renormalization group flow leads 
to spontaneous chiral symmetry breaking in vacuum, the symmetry gets restored in a second 



order phase transition at high temperature for vanishing quark mass. The description fn\ 
covers both the low temperature chiral perturbation theory domain of validity as well as the 
high temperature domain of critical phenomena. In particular, we obtain a precise estimate 



of the universal equation of state in the vicinity of critical points [n,18]. We explicitly 
connect the physics at zero temperature and realistic quark mass with the universal behavior 
near the critical temperature Tc and the chiral limit. An important property will be the 
observation that certain low energy properties are effectively independent of the details of the 
model even away from the phase transition. This behavior is caused by a strong attraction 
of the renormalization group flow to approximate partial infrared fixed points (infrared 
stability) [|I^,|I^. Within this approach at high density we find I^Dj a chiral symmetry 



restoring first order transition. As pointed out above these results imply the presence of a 
tricritical point with long-ranged correlations in the phase diagram. The lectures end with 
an outlook to densities above the chiral transition, where QCD is expected to behave as a 
color superconductor at low temperature. 

II. COARSE GRAINING AND NONPERTURBATIVE FLOW EQUATIONS 

A. Prom short to long distance scales 

Quantum chromodynamics describes qualitatively different physics at different length 
scales. The theory is asymptotically free [Q and at short distances or high energies the 
strong interaction dynamics of quarks and gluons can be determined from perturbation 
theory. On the other hand, at scales of a few hundred MeV confinement sets in and the 
spectrum of the theory consists only of color neutral states. The change of the strong gauge 
coupling as with scale has been convincingly demonstrated by a comparison of various 
measurements with the QCD prediction as shown in figure |l]. (Taken from [^.) The 



coupling increases with decreasing momenta Q and perturbation theory becomes invalid 
for Q ^ l.SGeV. Extrapolating QCD from short distance to long distance scales is clearly a 
nonperturbative problem. However, not only effective couplings but also the relevant degrees 
of freedom can change with scale. Low-energy degrees of freedom in strong interaction 



physics may comprise mesons, baryons and glueballs rather than quarks and gluons. Indeed, 
at low energies an essential part of strong interaction dynamics can be encoded in the masses 
and interactions of mesons. A prominent example of a systematic effective description of 
the low energy behavior of QCD is chiral perturbation theory^ [^]. It rather accurately 



describes the dynamics of the lightest hadronic bound states, i.e. the Goldstone bosons of 
spontaneous chiral symmetry breaking. 
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FIG. 1. Running of the strong gauge coupling by various types of measurements compared to 
theory. 

Spontaneous chiral symmetry breaking is an important building block for our under- 
standing of the hadron spectrum at low energies. The approximate chiral symmetry, hiding 
behind this phenomenon, originates in the fact that left- and right-handed components of 
free massless fermions do not communicate. In the limit of vanishing quark masses the 
classical or short distance QCD action does not couple left- and right-handed quarks. As 
a consequence it exhibits a global chiral invariance under [/^(A^/) x Uji{Nf) = SUL{Nf) x 
SUni^Nf) X Uy{1) X [/^(l). Here Nf denotes the number of massless quarks and the quark 
fields ip transform with different unitary transformations acting on the left- and right- 
handed components 
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^See also the lectures of this school from C.P. Burgess [23| for an introduction and references 



In reality, the quark masses differ from zero and the chiral symmetry is only an approximate 
symmetry. There are two especially light quark flavors with similar quark masses, rriu — rud- 
If the masses of u and d are taken to be the same they form an isospin doublet. The cor- 
responding approximate isospin symmetry of the QCD action is manifest at low energies in 
the observed pattern of bound states which occur in nearly degenerate multiplets: {p,n), 
{7c~^,TT^,7[~), . . . Even though the strange quark s is much heavier than u and d the asso- 
ciated symmetry for three degenerate flavors, termed the "eightfold way", can be observed 
as mesonic and baryonic levels grouped into multiplets of SU{3) — singlets, octets, decu- 
plets. The other flavors c, b, t remain singlets. If chiral symmetry was realized in the same 
manner, the energy levels would appear approximetely as multiplets of SUl(2) x SUr{2), 
or SUii?)) X SUb.{?>), respectively. The multiplets would then necessarily contain members 
of opposite parity which is not observed in the hadron spectrum. Furthermore the pion 
mass is small compared to the masses of all other hadrons. Together this indicates that 
the approximate chiral symmetry SUilNf) x SUptiNf) with Nf = 2 or 3 is spontaneously 
broken to the diagonal SUv{Nf) vector-like subgroup with the pions as the corresponding 
Goldstone bosons. 

SUL{Nf) X SUniNj) x Uv{l) x ^7^(1) 




spontaneously 
broken 



baryon 
number 
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chiral anomaly 



SUviNf) X Uvil) 

In addition, the axial abelian subgroup Ua{^) of the classical QCD action is broken in the 
quantum theory by an anomaly of the axial-vector current [^]. This breaking proceeds 
without the occurence of a Goldstone boson. The abelian Uv{l) subgroup corresponds to 
baryon number conservation. 



B. Kadanoff Wilson renormalization group 



A conceptually very appealing approach to bridge the gap between the short distance and 
the long distance physics relies on the general ideas of the Kadanoff- Wilson renormalization 
group 1 25, 26 1 . The renormalization group method consists in systematically reducing the 
number of degrees of freedom by integrating over short wavelength fluctuations. For an 
example, one may consider a spin system on a lattice with lattice spacing a. A possible 
strategy for integrating over short wavelength fluctuations is to form blocks of spins. In 
figure 2 the spins are grouped by fours and to each block one attributes an appropriate 
"average" spin. One may visualize such a procedure by imagining the spin system observed 
through a microscope with two different resolutions. At first one uses a relatively high 





FIG. 2. 



resolution, the mininiuin wavelength for fluctuations is ~ a. Subsequently one uses a longer 
wavelength, permitting the observation of details with dimension ^ 2a. The resolution 
becomes less since details having dimension between a and 2a have been integrated out. As 
a result one obtains an effective description for averages of spins on the sites of a coarse grain 
(block) lattice, incorporating the effects of the short wavelength fluctuations. One should 
note that for the averaging procedure (a -^ 2a) fluctuations on scales larger than the coarse 
graining length scale (2a) play no role. This observation can be used to obtain a relatively 
simple coarse graining description also in the case of a quantum fleld theory which will be 
considered in the following. 

C. Average action 



I will introduce here the average action r^ [^ which is based on the quantum fleld 
theoretical concept of the effective action F, i.e. the generating functional of the IPI Green 
functions. The fleld equations derived from the effective action include all quantum effects. 
In thermal and chemical equilibrium F includes in addition the thermal fluctuations and 
depends on the temperature and chemical potential. In statistical physics F corresponds to 
the free energy as a functional of some (space dependent) order parameter. 

• The average action F^ is a simple generalization of the effective action, with the dis- 
tinction that only fluctuations with momenta q^ ^ /c^ are included. 

This is achieved by implementing a momentum infrared cutoff ~ k into the functional 
integral which deflnes the effective action F. In the language of statistical physics, F^ is a 
coarse grained free energy with a coarse graining length scale ~ k~^. Lowering k results in a 
successive inclusion of fluctuations with momenta q^ ^ /c^ and therefore permits to explore 
the theory on larger and larger length scales. The average action F^ can be viewed as the 
effective action for averages of flelds over a volume with size fc~°' and is similar in spirit to 
the action for block-spins on the sites of a coarse lattice. By deflnition, the average action 
equals the standard effective action for fc = 0, i.e. Fq = F, since the infrared cutoff is absent 
and therefore all fluctuations are included. On the other hand, in a theory with a physical 
ultraviolet cutoff A we can associate Fa with the microscopic or classical action S since no 
fluctuations below A are effectively included. Thus the average action has the important 
property that 



• Tk interpolates between the classical action S and the effective action F as A; is lowered 
from the ultraviolet cutoff A to zero: lim Tu = S , lim Fi. = F . 

The ability to follow the evolution to A; ^ is equivalent to the ability to solve the theory. 
Most importantly, the dependence of the average action on the scale k is described by an 
exact nonperturbative flow equation which is presented in the next section. 

Let us consider the construction of F^ for a simple model with real scalar fields Xa, 
a = 1 . . . N, in d Euclidean dimensions with classical action S and sources J^. We start 
with the path integral representation of the generating functional for the connected Green 
functions 

W,[J] = InJDxexp i^-S[x] - ASk[x] + j d''xJa{x)x\x)^ , (2) 

where we have added to the classical action an infrared (IR) cutoff term AS'fc[x] which is 
quadratic in the fields and reads in momentum space 

^Su[x] = lJ-0y!M<l)Xa{-q)x%q). (3) 

Here the infrared cutoff function Rk is required to vanish for k ^ and to diverge for k -^ A 
and fixed g^. For A ^ oo this can be achieved, for example, by the exponential form 

For fiuctuations with small momenta g^ <^ A;^ this cutoff behaves as Rk{q^) ~ k'^ and allows 
for a simple interpretation. Since AS'fc[x] is quadratic in the fields, this means that all Fourier 
modes of x with momenta smaller than k acquire an effective mass ~ k. This additional mass 
term acts as an effective IR cutoff for the low momentum modes. In contrast, for q^ ^ A;^ the 
function Rk{q^) vanishes such that the functional integration of the high momentum modes 
is not disturbed. The term AS'fc[x] added to the classical action is the main ingredient for 
the construction of an effective action that contains all fluctuations with momenta q^ ^ k"^ 
whereas fluctuations with q^ ^ k"^ are suppressed. 

The expectation value of x, i-e. the macroscopic field 0, in the presence of AS'^fx] and J 
reads 

We note that the relation between and J is A;-dependent, (p = (pk{J) and therefore J = Jk{(f))- 
In terms of Wk the average action is deflned via a modifled Legendre transform 

Tfcb] = -Wk[J] + J d''xJa{x)r{x) - /\Su[<p\ (6) 

where we have subtracted the term AS'fc[0] on the r.h.s. This subtraction of the infrared 
cutoff term as a function of the macroscopic fleld is crucial for the definition of a reasonable 
coarse grained free energy with the property lim^^ArA: = S. It guarantees that the only 



-T— + (pRk 



difference between F^ and F is the effective infrared cutoff in the ffuctuations. 

Exercise: Check of the properties (i) hniFfc = F, (ii) hm F^ = S. 

(i) The first property follows immediately with linifc^o Rk = ^ from the absence of any IR 
cutoff term in the above (standard) definitions for the effective action. 

(ii) To establish the property Fa = 5* we consider an integral equation for F^ which is 
equivalent to (||). In an obvious matrix notation we use (^ 

exp {Wk[J]) = j Dx exp [-S[x] + j Jx-\jx Rkx) ■ 

Eliminating in this equation Wk and J with (^ 

-Wu[J]=Tk['i>] + ]^j <pRk<p- j J<P => J 

we obtain the integral equation 

exp(-FfcM) =jDxexp (-S[x] + / ^[x - 0])exp (-^[x - 0] Rk [x - <t>]) ■ (7) 

For k —* A the cutoff function R^ diverges and the term exp(— [x — (f)]Rk[x — </']/2) 
behaves as a delta functional ~ (5[x — 0], thus leading to the property F^ — > S* in this 
limit. 



Let us point out a few properties of the average action: 

1. All symmetries of the model which are respected by the IR cutoff term ASk are also 
symmetries of F^. In particular, this concerns translation and rotation invariance. 

2. The construction of the average action can be easily generalized to fermionic degrees 
of freedom. In particular, it is possible to incorporate chiral fermions since a chirally 
invariant cutoff Rk can be formulated p8 |. 



3. Gauge theories can be formulated along similar lines [P9|-|32| even though ASk may 



not be gauge invariant. In this case the usual Ward identities receive corrections for 



which one can derive closed expressions |^. These corrections vanish for k ^ 0. 



4. Despite the similar spirit one should note the difference in viewpoint to the Kadanoff- 



Wilson block-spin action p5| , P6| . The Wilsonian effective action realizes that physics 
with a given characteristic length scale / can be conveniently described by a functional 
integral with an ultraviolet cutoff A for the momenta where l~^ should be smaller than 
A, but not necessarily by a large factor. The Wilsonian effective action S"^ replaces 

10 



then the classical action in the functional integral. It is obtained by integrating out 
the fluctuations with momenta g^^A^. The n-point functions have to be computed 
from S"^ by further functional integration and are independent of A. 

In contrast, the average action F^ realizes the concept of a coarse grained free energy 
with a coarse graining length scale ~ fc~^. For each value of k the average action is 
related to the generating functional of a theory with a different action Sk = S + A^S. 
The n-point functions derived from F^ depend on k. To obtain /c-independent n- 
point functions their characteristic momentum scale l~^ should be much larger than 
k. The standard effective action F is obtained by following the flow for F^ to fc = 0, 
thus removing the infrared cutoff in the end. The Wilsonian effective action does not 



generate the IPI Green functions [^ 



D. Exact flow equation 

The dependence of the average action F^ on the coarse graining scale k is described by 



an exact nonperturbative flow equation [|1^] 

d 



di^'^ 




(2) I 




n'[<p] + Rk ^Rk\ . (8) 



Here t = ln(/c/A) denotes the logarithmic scale variable with some arbitrary momentum scale 
A. The trace involves only one integration as well as a summation over internal indices, and 
in momentum space it reads Tr = X^a / d'^q/i'^'^Y for the a = 1, . . . ,N component scalar fleld 
theory. The exact flow equation describes the scale dependence of F^ in terms of the inverse 

(2) 

average propagator F^ as given by the second functional derivative of F^ with respect to 
the fleld components 



Exercise: Derivation of the exact flow equation (|^). 
Let us write 

Ffc[0]=f,[0]-A5fc[0] (10) 

where according to (^ 

ffc[0] = -Wk[J]+Jd''xJ{x)<f){x) (11) 

and J = Jk{4>)- We consider for simplicity a one-component fleld and derive flrst the 
scale dependence of F: 

11 



lf,,H-(^)l,-/.^.i;l^./.^.„.)ME). ,, 



With (f){x) = 6Wk/SJ{x) the last two terms in (|T2]) cancel. The t-derivative of Wk is 
obtained from its defining functional integral (|]) and yields 

-f.M = (^A5,M) = {- j d'xdSx{x)j^Mx.y)x{y)) ■ (13) 

where Rk{x,y) = Rk{—d1)5{x — y). Let G{x,y) = 6^Wk/6J{x)6J{y) denote the con- 
nected 2-point function and decompose 

{X{x)x{y)) = G{x, y) + (x(x))(x(y)) = G{x, y) + <P{x)<P{y) . (14) 

Plugging this decomposition into (pTsD the scale dependence of F^ can be expressed as 

-gflA't'] = 2 / d^'^'i'v \ G(x,y)igRk(x, y) + il>(x)—Rt(x,y)4,{y) I 
The exact flow equation for the average action F^ follows now with ( p!OD 



^FJ01 = -Tr|G^i?fc 
dt ^^^^ 2 1 9t 



in. 

2 



^(2)[ 



-1 d 



n'[^] + Rk\ ^^Rk\ (16) 



where we have used that F^ {x,y) = 5'^Tk/S(j){x)6(j){y) = 6J{x)/6(f){y) is the inverse of 
G{x,y) = 6'^Wk/SJ{x)6J{y) = 6(t>{x)/6J{y) to obtain the last equation. 



Let us point out a few properties of the exact flow equation: 

/"ON 

1. The flow equation (|) closely resembles a one-loop equation. Replacing F^ by the 
second functional derivative of the classical action, S^^\ one obtains the corresponding 
one-loop result. Indeed, the one-loop formula for F^ reads 

Ffc[0] = 5H + ^ Tr In (S^^) [0] + R^) (17) 



and taking a t-derivative of ([17D gives a one-loop flow equation very similar to (^. 
The "renormalization group improvement" S*^^^ — > F^ turns the one-loop flow equa- 
tion into an exact nonperturbative flow equation which includes the effects from all 

(2) 

loops. Replacing the propagator and vertices appearing in F^ by the ones derived 
from the classical action, but with running fc-dependent couplings, and expanding 
the result to lowest non-trivial order in the coupling constants one recovers standard 
renormalization group improved one-loop perturbation theory. 
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2. The additional cutoff function Rk with a form hke the one given in eq. (Q) renders the 
momentum integration imphed in the trace of (§) both infrared and ultraviolet finite. 
In particular, the direct implementation of the additional mass-like term R^ ~ k"^ 
for q^ <^ k'^ into the inverse average propagator makes the formulation suitable for 
dealing with theories which are plagued by infrared problems in perturbation theory. 
We note that the derivation of the exact flow equation does not depend on the par- 
ticular choice of the cutoff function. Ultraviolet finiteness, however, is related to a 
fast decay of dtRk for q^ S> k'^. If for some other choice of Rk the right hand side of 
the flow equation would not remain ultraviolet finite this would indicate that the high 
momentum modes have not yet been integrated out completely in the computation of 
Ffc. Unless stated otherwise we will always assume a sufficiently fast decaying choice 
of Rk in the following. 

Of course, the particular choice for the infrared cutoff function should have no effect 
on the physical results for /c — >■ 0. Different choices of Rk correspond to different 
trajectories in the space of effective actions along which the unique infrared limit Fq 
is reached. Nevertheless, once approximations are applied not only the trajectory 
but also its end point may depend on the precise definition of the function Rk- This 
dependence may be used to study the robustness of the approximation. 



3. Flow equations for n-point functions can be easily obtained from (^ by differentiation. 
The flow equation for the two-point function F^ involves the three and four-point 
functions, F^ and F^ \ respectively: 



^rf(.,.) 



dt d(j){—q)d(f){q) 
dRk d 



dt d(f){—q) 




(-1) 



+ Rk 



F 



^(3) 



(2) 



Rk 



J- f. 



Fl^) + Rk 



Tu + Rk 



-1 



^(3) 



+ Rk 



^(4) 



Ff + Rk 



-1 



Fr. + Rk 



-1 



(18) 



(n) 



.(n+1) 



In general, the flow equation for F^ involves F^ and F 



■^in+2) 



We emphasize that the flow equation (^ is equivalent to the Wilsonian exact renor- 
malization group equation |26|j3^-p8| . The latter describes how the Wilsonian effective 



action S*^ changes with the ultraviolet cutoff A. Polchinski's continuum version of the 



Wilsonian flow equation |^ can be transformed into eq. 
transform and a suitable variable redefinition pU| , ^ .p| 



by means of a Legendre 



5. Extensions of the flow equations to gauge fields 12 
able. 



2i and fermions 28,42 are avail 



^See also the contribution from C. Kim [Elf| this school 
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E. Elements of effective field theory 

A strict derivation of an effective low energy description from QCD is still missing. 
Yet, predictions from effective low energy models often show convincing agreement with the 
results of real and numerical experiments. Let us consider some important aspects for the 
success of an effective field theory description. They find a natural theoretical basis in the 
framework of the average action. 

1) Decoupling of "heavy" degrees of freedom — At any scale k only fluctuations 
with momenta in a small range around k influence the renormalization group flow of the 
average action F^. This expresses the fact that the momentum integration implied by 
the trace on the r.h.s. of the exact flow equation (H) is dominated by momenta g^ ~ fc^, 
schematically 



or, 1 



dk 2 





with the full A;-dependent propagator associated to the propagator line and the dot denotes 
the insertion dtRk- The effects of the high momentum modes with q'^^k'^ determine the 
precise form of the average action Tk at the scale k. Modes with g^ ^ k"^ have been "integrated 
out". In particular, for any given k = ko their only effect is to determine the initial value 
Ffcg for the solution of the flow equation for the low momentum fluctuations with q"^ ^k^. 

Most importantly, for each range of k only those degrees of freedom have to be included 
which are relevant in the corresponding momentum range. In QCD these may comprise 
compounds of quarks and gluons. "Heavy" degrees of freedom effectively decouple from the 
flow of Ffc once k drops below their mass m, since m represents a physical infrared cutoff for 
fluctuations with momenta g^ ~ k"^ ^rn^ . We will observe in section piBj the occurence of 
mass threshold functions, which explicitly describe the decoupling of heavy modes, as one 
of the important nonperturbative ingredients of the flow equations. 

As a prominent example for an effective field theory, chiral perturbation theory describes 
the IR behavior of QCD in terms of the lightest mesons, i.e. the Goldstone bosons of spon- 
taneous chiral symmetry breaking. This yields a very successful effective formulation of 
strong interactions dynamics for momentum scales up to a few hundred MeV. For some- 
what higher scales additional degrees of freedom like the sigma meson or the light quark 
flavors will become important and should be included explicitly. The linear quark meson 
model based on these degrees of freedom will be introduced in the next section. 

2) Infrared stability — The predictive power of an effective low energy description 
crucially depends on how sensitively the infrared value F = linifc^o Tfc depends on the initial 
value Ffcg. Here fco plays the role of an ultraviolet cutoff scale of the low energy description. 
Indeed, as the coarse graining scale k is lowered from fco to zero, the "resolution" is smeared 
out and the detailed information of the short distance physics can be lost. (On the other 
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hand, the "observable volume" is increased and long distance aspects such as collective 
phenomena become visible.) 

There is a prominent example for insensitivity of long distance properties to details at 
short distances: Systems of statistical mechanics near critical points, where they undergo a 
second order phase transition. For the example of a spin system near its critical temperature 
the equation of state, which relates for a given temperature the magnetization to an external 
magnetic field, is independent of the microscopic details up to the short distance value of 
two parameters. These can be related to the deviation from the critical temperature and 
to the deviation from a nonzero magnetic field. Stated differently, only two parameters of 
the short distance effective action F^q have to be "finetuned" in order to be at the phase 
transition. These few relevant parameters are typically accompanied by many irrelevant 
parameters which parametrize F^.Q As the coarse graining scale k is lowered from ko to zero, 
the running (appropriately rescaled) irrelevant parameters are attracted to a fixed point, 
whereas relevant parameters are driven away from this point. Figure 3 shows, schematically, 
the vicinity of a fixed point for the case of one irrelevant parameter and one relevant param- 
eter. The arrows indicate the renormalization group flow. Since the irrelevant parameter is 
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FIG. 3. 

driven to the same fixed point irrespective of its initial value (once the relevant parameter is 
tuned to criticality), the short distance information about it gets lost in the infrared. This 
behavior is crucial for the phenomenon of universality for critical phenomena (cf. section 

rvTi). 



We will observe within the linear quark meson an important different example of the 
insensitivity of physics at long distances to short distance details even away from the phase 
transition. A crucial observation is the strong attraction of the renormalization group flow 



to approximate fixed points. We will discuss this remarkable property in section |TTB. 

3) Symmetries — Symmetries constrain the possible form of an effective model. It 
is an important property that the coarse grained free energy F^ respects all symmetries of 



^For a more detailed introduction and a systematic classification of operators as relevant, marginal 
or irrelevant see e.g. ref. 0]. 
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the model, in particular rotation and translation symmetries. In consequence, Tk can be 
expanded in terms of invariants with respect to the symmetries with couplings depending on 
k. We will frequently employ an expansion of r^ where the invariants are ordered according 
to the number of derivatives of the fields. For the example of the 0(iV)-symmetric scalar 
theory this yields (a = 1, . . . , A^) 

Tfc = Jd''x |f/,(p) + ^Zk{p)d'^<j)ad,r + \Yk{p)d>'pd^p + 0{d^)^ (19) 

with p = 4>"'(f)a- Here Uk{p) corresponds to the most general 0(A^)-symmetric non-derivative, 
potential term. The derivative terms contain a field dependent wave function renormalization 
factor Zk{p) plus a function Vfc(p) accounting for a possible different index structure of the 
kinetic term for N > 2. Going further would require the consideration of terms with four 
derivatives and so on. The (approximate) chiral symmetry of QCD will play an important 
role for the construction of the quark meson model which will be discussed in the following. 



III. THE LINEAR QUARK MESON MODEL 

Before discussing the nonzero temperature and density behavior of strong interaction 
physics we will review some of its zero temperature features. This will be done within the 
framework of a linear quark meson model as an effective description for QCD for scales below 
the mesonic compositeness scale of approximately A;$ ^ 600 MeV. Relating this model to 
QCD in a semi-quantitative way in subsection |IIIA| will allow us to gain some information 
on the initial value for the effective average action at the scale A;$. We emphasize, however, 
that the quantitative aspects of the derivation of the effective quark meson model from QCD 
will not be relevant for our practical calculations in the mesonic sector. This is related to 
the infrared stability mentioned in the previous section and which will be made quantitative 



in IIIB. 



A. A short scale history 

We give here a brief semi-quantitative picture of the relevant scales and the physical 
degrees of freedom which appear in relation to the phenomenon of chiral symmetry breaking 
in vacuum. See also the reviews in P^JT7|] . Most of this will be explained in more detail 
below whereas other parts are rather well established features of strong interaction physics. 

We will distinguish five qualitatively different ranges of scales: 

1. At sufficiently high momentum scales, say, 

k^kp- 1.5 GeV 

the relevant degrees of freedom of strong interactions are quarks and gluons and their 
dynamics is well described by perturbative QCD. 
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2. For decreasing momentum scales in the range 

k^ ~ 600 MeV ^ A; ^ A;^ ~ 1.5 GeV 

the dynamical degrees of freedom are still quarks and gluons. Yet, as k is lowered part 
of their dynamics becomes dominated by effective non-local quark interactions which 
cannot be fully accessed perturbatively. 

3. At still lower scales this situation changes dramatically. Quarks and gluons are sup- 
plemented by mesonic bound states as additional degrees of freedom which are formed 
at a scale /c$ ~ 600 MeV. We emphasize that /c$ is well separated from Aqcd — 
200 MeV where confinement sets in and from the constituent masses of the quarks 
Mq ~ (300 — 350) MeV. This implies that below the compositeness scale A;$ there 
exists a hybrid description in term of quarks and mesons. It is important to note 
that for scales not too much smaller than A;$ chiral symmetry remains unbroken. This 
situation holds down to a scale k^sB — 400 MeV at which the scalar meson potential 
develops a non-trivial minimum thus breaking chiral symmetry spontaneously. The 
meson dynamics within the range 

k^SB ^ 400 MeV ^ A; ^ A;$ ^ 600 MeV 

is dominated by light current quarks with a strong Yukawa coupling /i^/(47r) ^ as{k) 
to mesons. We will thus assume that the leading gluon effects are included below A;$ 
already in the formation of mesons. Near k^sB also fluctuations of the light scalar 
mesons become important as their initially large renormalized mass approaches zero. 
Other hadronic bound states like vector mesons or baryons should have masses larger 
than those of the lightest scalar mesons, in particular near k^sB, and give therefore 
only subleading contributions to the dynamics. This leads us to a simple linear model 
of quarks and scalar mesons as an effective description of QCD for scales below kq>. 

4. As one evolves to scales below k^sB the Yukawa coupling decreases whereas ag in- 
creases. Of course, getting closer to Aqcd it is no longer justified to neglect in the 
quark sector the QCD effects which go beyond the dynamics of the simple effective 
quark meson model. On the other hand, the final IR value of the Yukawa coupling 
h is fixed by the typical values of constituent quark masses Mq ~ 300 MeV to be 
/i^/(47r) ~ 3.4. One may therefore speculate that the domination of the Yukawa inter- 
action persists even for the interval 

Mq ~ 300 MeV ^k^ k^sB ^ 400 MeV 

below which the quarks decouple from the evolution of the mesonic degrees of freedom 
altogether. Of course, details of the gluonic interactions are expected to be crucial 
for an understanding of quark and gluon confinement. Strong interaction effects may 
dramatically change the momentum dependence of the quark propagator for k and g^ 
around Aqcd- Yet, there is no coupling of the gluons to the color neutral mesons. As 
long as one is only interested in the dynamics of the mesons one is led to expect that 
confinement effects are quantitatively not too important. 
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5. Because of the effective decoupling of the quarks and therefore of the whole colored 
sector the details of confinement have only little influence on the mesonic dynamics 
for scales 

k^MgC^ 300 MeV . 

Here quarks and gluons disappear effectively from the spectrum and one is left with 
the pions. For scales below the pion mass the flow of the couplings stops. 

1. 'Integrating out gluons' 

Let us now discuss the above picture with different ranges of scales in more detail. In 
order to obtain the effective action at the compositeness scale fc$ from short distance QCD 
two steps have to be carried out. In the flrst step one computes at the scale kp ~ 1.5 GeV an 
effective action involving only quarks. This step integrates out the gluon degrees of freedom 
in a "quenched approximation". More precisely, one solves a truncated flow equation for 
QCD with quark and gluon degrees of freedom in presence of an effective infrared cutoff 
kp ^ 1.5 GeV in the quark propagators. This procedure is outlined in |H5|. The exact flow 



equation to be used for this purpose is obtained by lowering the infrared cutoff Rk for the 
gluons while keeping the one for the quarks flxed. Subsequently, the gluons are eliminated 
by solving the fleld equations for the gluon flelds as functionals of the quarks. This will 
result in a non-trivial momentum dependence of the quark propagator and effective non- 
local four and higher quark interactions. Because of the infrared cutoff kp the resulting 
effective action for the quarks resembles closely the one for heavy quarks. The dominant 
effect is the appearance of an effective quark potential (similar to the one for the charm 
quark) which describes the effective four-quark interactions. First promising results within 
this approach include an estimate of the heavy quark effective potential valid for momenta 
y/q^^300 — 500 MeV |4^,47]. The inverse quark propagator is found in this computation to 



remain very well approximated by the simple classical momentum dependence 'y/j.q^. 

2. Formation of mesonic bound states 

In the second step one has to lower the infrared cutoff in the effective non-local quark 
model in order to extrapolate from kp to fc<j,. This task can be carried out by means of 
the flow equation for quarks only, starting at kp with an initial value F^p ['?/'] as obtained 
after integrating out the gluons. For fermions the trace in (|^) has to be replaced by a 
supertrace in order to account for the minus sign related to Grassmann variables [^. A 



flrst investigation in this direction ^^ has used a chirally invariant four quark interaction 



'dressed" one-gluon exchange) whose most general momentum dependence was retained 




' S{pi + P2 - P3 - P4))^i''\pi,P2,P3,PA)M{pi,P2,P3,P4) , (20) 

M{p„P2,P3,P4) = - {^P:{-p^)Y{n,'^P]{-p3)} {^',{p,h,{T.),Ve{P2)} ■ (21) 



The curled brackets indicate contractions over spinor indices, i,j, k, I = l...Nc are the colour 

indices and a,b = l.-.A^/ the flavour indices of the quarks. By an appropriate Fierz trans- 
formation and using the identity 

iTnAT,)k=ls-Si-^fA (22) 
one can split Ai into three terms 

M = M^ + Mp + Mp (23) 

\ {'A:(-Pi)7'^- (P2)} {^^■(P4)7'^J(-P3)} (24) 
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M, = -^ {v5:(-pi)z7mC(-P3)} {i^liP4)t7'^',iP2)} ■ (26) 

In terms of the Lorentz invariants 

S = (pi +P2f = {Pi+Pif 

t={pi~P^f = {p2-PAf (27) 

we recognize that the quantum numbers of the fermion bilinears in A^o- correspond to colour 
singlet, flavour non-singlet scalars in the s-channel and similarly for spin-one mesons for M.p. 



Following [^] we associate these terms with the scalar mesons of the linear a-model and with 
the p-mesons. The bilinears in the last term M.p correspond to a colour and flavour singlet 
spin-one boson in the t-channel. These are the quantum numbers of the pomeron. In the 
following we neglect interactions in the vector meson and pomeron channels and only retain 
the contribution M.^- We will discuss this approximation in more detail in section p.V D| . 
The matrices m and ifi are hermitian and m + irfi^^ forms therefore the most general quark 
mass matrix. (Our chiral conventions |^ where the hermitean part of the mass matrix 
is multiplied by 75 may be somewhat unusual but they are quite convenient for Euclidean 
calculations.) With V{q^) the heavy quark potential in a Fourier representation, the initial 
value at fcp = 1.5 GeV was taken as {Z^^^ = Z^^^ip"^ = ~kp)) 

>.it\p.P.,P.P.)Zil = iv((p. - P,r) = (^ + i^^^^ . (28) 

This corresponds to an approximation by a one gluon exchange term ~ as{kp) and a string 
tension A ~ 0.18 GeV^ and is in reasonable agreement with the form computed recently EH 



from the solution of flow equations. In the simplifled ansatz (p8|) the string tension introduces 
a second scale in addition to kp and indeed the incorporation of gluon fluctuations is a crucial 
ingredient for the emergence of mesonic bound states. For a more precise treatment |^ of 
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the four-quark interaction at the scale fc$ this second scale is set by the running of a^ or 

^QCD- 

The evolution equation for the function A^ for k < kp can be derived from the fermionic 
version of (§) and the truncation (0). Since A^ depends on six independent momentum 
invariants it is a partial differential equation for a function depending on seven variables 
and has to be solved numerically |^ . The "initial value" (pS]) corresponds to the t-channel 
exchange of a "dressed" colored gluonic state and it is exciting to realize that the evolution of 
A^ leads at lower scales to a momentum dependence representing the exchange of colorless 
mesonic bound states. At the compositeness scale 

k^ ~ 630 MeV (29) 

one finds [^ an approximate factorization 

Ait^(Pi,P2,P3,P4) = 9{PuP2)G{s)g{p3,p4) + ... (30) 

which indicates the formation of mesonic bound states. Here g{pi,P2) denotes the amputated 
Bethe-Salpeter wave function and G{s) is the mesonic bound state propagator displaying 
a pole-like structure in the s-channel if it is continued to negative s = (pi + ^2)^- The 
dots indicate the part of A^ which does not factorize and which will be neglected in the 
following. In the limit where the momentum dependence of g and G is neglected we recover 
the four-quark interaction of the Nambu-Jona-Lasinio model [^S|-pT|. 



3. Quark meson model 

For scales below the mesonic compositeness scale kij, a description of strong interaction 
physics in terms of quark fields alone would be rather inefficient. Finding physically rea- 
sonable truncations of the effective average action should be much easier once composite 
fields for the mesons are introduced. The exact renormalization group equation can indeed 
be supplemented by an exact formalism for the introduction of composite field variables or, 



more generally, a change of variables [^. For our purpose, this amounts in practice to 
inserting at the scale fc$ the identities 

1 = const / VffA 



X exp j- tr (ai - K\G -m\- O^G) i (a a - GKa -tha- GO) \ 

1 = const / Vau (31) 



X exp |- tr Uh - KliG - ni^n - O^^^^G) i Uh - GKh - m^ - GC^^^ 

into the functional integral which formally defines the quark effective average action. Here we 
have used the shorthand notation A^GB = J j^;4aAl{q)G"'^{q)Bb{q), and Ka,h are sources 
for the collective fields o'a,h which correspond in turn to the anti-hermitian and hermitian 
parts of the meson field $. They are associated to the fermion bilinear operators 0[iIj], 
O^^^lip] whose Fourier components read 
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o^'\iq) = -j -0j-M-p^p + <iWii>)i^Mp + q) ■ (32) 

The choice of g{—p,p + q) as the bound state wave function renormahzation and of G{q) 
as its propagator guarantees that the four-quark interaction contained in (|3TD cancels the 
dominant factorizing part of the QCD-induced non-local four-quark interaction Eqs.(pl|), 
(0). In addition, one may choose 

^lab = m,,(0)(7-i(0,0)Z^,fc,(0) (33) 

such that the explicit quark mass term cancels out for g = 0. The remaining quark bilinear 
is ~ m{q) — m{0)Z^^k^{0)g{—q, q)/[Z^^k^{q)g{0, 0)]. It vanishes for zero momentum and will 
be neglected in the following. Without loss of generality we can take m real and diagonal 
and m = 0. 

In consequence, we have replaced at the scale fc$ the effective quark action (|21|) with 
(0) by an effective quark meson action given by 

ffc = rfc-i|dStr($tj + jt$) 

Tk = Jd^xUk{<^,<^^) (34) 

"^ ^ ' Z$,,(g)g2tr l^^q^q)] + Z^MMqh'q.riq) 



(27r)'^ k 



(27r)' 
X r{q) (^^$a.(p) - ^^$L(-P)) ^\q-p) 



At the scale fc$ the inverse scalar propagator is related to G{q) in (pOD by 

G'\q') = 2ml+2Z^,k,iq)q'- (35) 

This fixes the term in Uk^ which is quadratic in $ to be positive, Uk^ = rn\ tr$''"$ + . . .. 
The higher order terms in Uk^ cannot be determined in the approximation (|2T|) since they 
correspond to terms involving six or more quark fields. The initial value of the Yukawa 
coupling corresponds to the "quark wave function in the meson" in (^), i.e. 

hk^{-q,q-p) = g{-q,q-p) (36) 

which can be normalized with hk^{0,0) = 5'(0,0) = 1. We observe that the explicit chiral 
symmetry breaking from non-vanishing current quark masses appears now in the form of a 
meson source term with 

J = 2ml^Z^^k^{0)g'^{0, 0) (mab + irfiab) = 2Z^^k^ml^dieig{mu, md, rus) . (37) 
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This induces a non- vanishing ($) and an effective quark mass Mq through the Yukawa 
couphng. We note that the current quark mass nig and the constituent quark mass Mg ~ 
hk ($) are identical at the scale /c$. (By solving the field equation for <l> as a functional of 
ip, ip (with Uk = m|tr$^$) one recovers from (|3^) the effective quark action (^ll) . For a 
generalization beyond the approximation of a four-quark interaction or a quadratic potential 



see ref ^j.) Spontaneous chiral symmetry breaking can be described in this language by a 
non- vanishing (<l>) in the limit j —>■ 0. 

The effective potential Uk{^) must be invariant under the chiral SUl{N) x SUii{N) flavor 
symmetry. In fact, the axial anomaly of QCD breaks the Abelian Ua_{^) symmetry. The 
resulting Ua{^) violating multi-quark interactions^ lead to corresponding Ua{^) violating 
terms in Uk{^)- Accordingly, the most general effective potential Uk is a function of the 
A^ + 1 independent C and P conserving SUl{N) x SUr{N) invariants 

p = tr^"^*!) , 

r, ~ trf^t^-lpV , t = 2,...,N , (38) 



A^ 
^ = det $ + det <l>^ . 

We will concentrate in this work on the two flavor case (A^ = 2) and comment on the 



effects of including the strange quark in section [IV D| . Furthermore we will neglect isospin 
violation and therefore consider a singlet source term j proportional to the average light 
current quark mass m = ^{niu + md)- Due to the f/^(l)-anomaly there is a mass split for 
the mesons described by $ 

<^ = \{a-ir^') + ]^(j' + i7i^)Tk. (39) 

The scalar triplet (oq) and the pseudoscalar singlet {j]') receive a large mass whereas the 
pseudoscalar triplet (vr) and the scalar singlet (a) remains light. From the measured values 
rrin' , m^Q ~ 1 GeV it is evident that a decoupling of these mesons is presumably a very 
realistic limit. (In thermal equilibrium at high temperature this decoupling is not obvious. 



We will comment on this point in section |I V D| . ) It can be achieved in a chirally invariant 



way and leads to the well known 0(4) symmetric Gell-Mann-Levy linear sigma model ||5^ 
which is, however, coupled to quarks now. This is the two flavor linear quark meson model 
which we will study in the next sections. For this model the effective potential Uk is a 
function of p only. 

The quantities which are directly connected to chiral symmetry breaking depend on 
the fc-dependent expectation value ($)^ = a^^k as given by the minimum of the effective 
potential 



^A first investigation for the computation of tfie anomaly term in the fermionic average action 



can be found in |53]. 
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In terms of the renorinalized expectation value 

o"o,fe = Z^f^'ao^k (41) 

we obtain the following expressions for quantities as the pion decay constant /^, chiral con- 
densate (ipip), constituent quark mass Mg and pion and sigma mass, m^r and mo-, respectively 
(d = 4) 



f'lT,k 


= 2(To,fc , 






~^^)k 


= -2m^, 


' ry-l/2 


) 


M,,k 


= hkCo^k , 


ml,k 


7-l/2^L> 7-1/2 J 

Cro,fc '^<yo,k 


Kk 


_ ^-1/2^ 


^*"+4A,ao% 





(42) 



Here we have defined the dimensionless, renormalized couplings 



dp^^'--'''^^ (43) 



Afc — ^$,fc n_9 (P ~ ^CTq I,; , 



/it: — -^^rts I, -^!^,;, i-Zi 



We are interested in the "physical values" of the quantities (|^) in the limit k -^ where 

V,fc=0' 



the infrared cutoff is removed, i.e. /^r = fn,k=o, fn'i = fn^ k=oy ^^'^■ 



4- Initial conditions 

At the scale kip the propagator G and the wave function g{—q, q — p) should be optimized 
for a most complete elimination of terms quartic in the quark fields. In the present context 
we will, however, neglect the momentum dependence of Z^^^, Zq,^^ and h^. We will choose a 
normalization of ip, $ such that Z^^^^ = h^^ = 1. We therefore need as initial values at the 
scale fc# the scalar wave function renormalization Z^ ^^ and the shape of the potential Uk^ ■ 
We will make here the important assumption that Zq,^k is small at the compositeness scale 
/c$ (similarly to what is usually assumed in Nambu-Jona-Lasinio-like models) 

Z^,k^ < 1 . (44) 

1/9 1 — 

This results in a large value of the renormalized Yukawa coupling hk = Z^J Z^f^hk. A 
large value of h^^ is phenomenologically suggested by the comparably large value of the 
constituent quark mass Mq. The latter is related to the value of the Yukawa couphng for 
/c — s- and the pion decay constant /^ = 92.4 MeV by Mq = hf.„/2 (with h = hk=o), and 
Mq ~ 300 MeV implies h'^/Air ~ 3.4. For increasing k the value of the Yukawa coupling grows 
rapidly for k^Mq. Our assumption of a large initial value for hk^ is therefore equivalent 
to the assumption that the truncation (0) can be used up to the vicinity of the Landau 
pole of hk- The existence of a strong Yukawa coupling enhances the predictive power of our 
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approach considerably. It implies a fast approach of the running couplings to partial infrared 
fixed points as is shown in section |1IIB [19,17|. In consequence, the detailed form of Uk^ 



becomes unimportant, except for the value of one relevant parameter corresponding to the 
scalar mass term m|^. In this work we fix m|_^ such that /^ = 92.4 MeV for m^ = 135 MeV. 
The value /,r = 92.4 MeV (for 171^^ = 135 MeV) sets our unit of mass for two flavor QCD 
which is, of course, not directly accessible by observation. Besides m^^^ (or /tt) the other 
input parameter used in this work is the constituent quark mass Mg which determines the 
scale fc$ at which h^^ becomes very large. We consider a range 300 MeV^Mg^350 MeV and 
find a rather weak dependence of our results on the precise value of Mg. 

All quantities in our truncation of F^ are now fixed and we may follow the fiow of r^ 
to k -^ 0. In this context it is important that the formalism for composite fields |Q also 



provides an infrared cutoff in the meson propagator. The flow equations are therefore exactly 
of the form (|^), with quarks and mesons treated on an equal footing. At the compositeness 
scale the quadratic term of Uk^ = rn\^ Tr $^$ + . . . is positive and the minimum of Uk^ 
therefore occurs for $ = 0. Spontaneous chiral symmetry breaking is described by a non- 
vanishing expectation value ($) in absence of quark masses. This follows from the change 
of the shape of the effective potential Uk as k flows from /c$ to zero. The large renormalized 
Yukawa coupling rapidly drives the scalar mass term to negative values and leads to a 
potential minimum away from the origin at some scale fc^ge < ^* such that flnally ($) = 



ctq 7^ for A; ^ [0,|T^,|l^]. This concludes our overview of the general features of chiral 



symmetry breaking in the context of flow equations. 

B. Flow equations and infrared stability 

1. Flow equation for the effective potential 

The dependence of the effective action Tk on the infrared cutoff scale k is given by an 
exact flow equation (Pf), which for fermionic flelds if) (quarks) and bosonic flelds $ (mesons) 



reads [|T6],|8l {t = \n{k/k^)) 



|r.|^.,*l = Ixr {^ (rf |V,.*1 + «.)-) - Tr {^ (r?'|^.,*| H- «,)-'} . (45) 
Here F^ is the matrix of second functional derivatives of F^ with respect to both fermionic 



and bosonic fleld components. The flrst trace on the right hand side of ( [45D effectively 
runs only over the bosonic degrees of freedom. It implies a momentum integration and a 
summation over flavor indices. The second trace runs over the fermionic degrees of freedom 
and contains in addition a summation over Dirac and color indices. The infrared cutoff 
function Rk has a block substructure with entries RkB and RkF for the bosonic and the 
fermionic flelds, respectively. 

We compute the flow equation for the effective potential Uk from equation (|45|) using the 
ansatz (BSI) for F^. The flow equation has a bosonic and fermionic contribution 



§-^Uk{p) = ^^Uksip) + ^^UkFip) . (46) 
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Let us first concentrate on the bosonic part by neglecting for a moment the quarks and 
compute dUks/dt . The effective potential Uk is obtained from F^ evaluated for constant 
fields. The bosonic contribution follows as 

In . ^ ^ f dS dR,Biq') j 3 1 1 

dt '""^P' 2 J (27r)4 dt \ Z^^,P,B{q^) + f/^(p) ^ Z^^.P.siq^) + U'^{p) + 2pt/^(p) j ' ^ ^ 

Here primes denote derivatives with respect to p and one observes the appearance of the 
(massless) inverse average propagator 

PkB{q^) = q^ + Z^^Rkiq') = ^ _ f^^y,, • (48) 

For p different from zero one recognizes the first term in (^) as the contribution from the 
pions, the Goldstone bosons of chiral symmetry breaking. (The mass term f/^ vanishes at the 
minimum of the potential.) The second contribution is then related to the radial or a-mode. 



Exercise: To obtain dUks/dt we calculate the trace in (^) for small field fiuctuations 
around a constant background configuration, which we take without loss of generality to be 
~ (t)5ai and p = 0^/2. The inverse propagator F^ can be computed from our ansatz ( |3l[ ) 
by expanding U^ to second order in small fiuctuations x, ^a{x) = (j)6ai + Xa{x) 

u, = u,{p) + uUp) (^l<i>a{x)<i>%x) - p) + ^uM (^^M^Wi^) - p)' + • • • 

= U,{p) + U',{p) (^XaX'^ + 0Xi) + lu'M<P\l + ... (49) 

and we obtain 

(rl'^)^^ (g, q') = {Z,q' + Uiip) + 2pU'M) {2^)%q - q') 

which yields (^T]). 



For the study of phase transitions it is convenient to work with rescaled, dimensionless 
and renormalized variables. We introduce 

u{trp) = k-''Uk{p), p = Z^^,k'~''p, hk = Zl]l'Z^lk''-^hk. (51) 

With 

-u{t, p)|p = -du{t, p) + {d-2 + r]^)pu'{t, p) + k-'^-U (p(p))|^ (52) 

one obtains from ( ^71) the evolution equation for the dimensionless potential. Here the 
anomalous dimension 77$ arises from the t-derivative acting on Zk and is given by 
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V^ = -yT(ln^#,fc) . (53) 

The fermionic contribution to the evolution equation for the effective potential can be com- 
puted without additional effort from the ansatz (|34D since the fermionic fields appear only 
quadratically. The respective flow equations is obtained by taking the second functional 
derivative evaluated at ip = 4' = 0- Combining the bosonic and the fermionic contributions 



one obtains the flow equation |T7 

d 

-—u = —du + (c? — 2 + ?7$) pu' 

9t (54) 



+ 2v, {3/o'(u ; v^) + I'oiu' + 2pu"; r/$) - 2-2^'N4^^\^ph'; 



W. 



Here v^^ = 2'^+^7r^/^r((i/2) is a prefactor depending on the dimension d and primes now 
denote derivatives with respect to p. The first two terms of the second line in ( ^4|) denote the 
contributions from the pions and the a-resonance (cf. (0)) and the last term corresponds 
to the fermionic contribution from the u, d quarks. The number of quark colors will always 
be A'^c = 3 in the following. 

The symbols /^, l^f'''^ in (Q) denote bosonic and fermionic mass threshold functions which 
contain the momentum integral implied in the traces on the r.h.s. of the exact flow equation 
(^). The threshold functions describe the decoupling of massive modes and provide an 
important non-perturbative ingredient. For instance, the bosonic threshold functions read 

l<i(yj. r,^) - ^ + ^n,o ij2n-d f d'^q 1 dRk 1 . . 

"^'''^*^" 4 ""'^ i (2vr)^Z^,. 9t [P,^(g2) + p^]«H-i- ^^^) 

These functions decrease ~ ti)^("+i) for w ':$> 1. Since typically w = M"^ /k"^ with M a mass, 
the main effect of the threshold functions is to cut off fluctuations of particles with masses 
M^ ^ k"^. Once the scale k is changed below a certain mass threshold, the corresponding 
particle no longer contributes to the evolution and decouples smoothly. 

Eq. (^) is a partial differential equation for the effective potential u{t,p) which has to 
be supplemented by the flow equation for the Yukawa coupling hk and expressions for the 
anomalous dimensions, where 

rj^ = -T:(}^^^,k) ■ (56) 

The corresponding flow equations and further details can be found in ref. [|17|,|19|. We will 
consider them in the next section in a limit where they can be solved analytically. We note 
that the running dimensionless renormalized expectation value k = 2A;^"'^Z$ ^aQ^,, with Wo^k 
the fc-dependent expectation value of $, may be computed for each k directly from the 
condition 



u'it,K) = ^k-'^Z^^J\ (57) 

\/2k, 
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2. Infrared stability 



Most importantly, one finds that the system of flow equations for the effective potential 
Uk{p), the Yukawa couphng hk and the wave function renormahzations Z^^^, Z^^^ exhibits 



an approximate partial fixed point [|T^,|l^. The small initial value (|4^) of the scalar wave 
function renormalization Z$ ^^ at the scale fc$ results in a large renormalized meson mass 
term Z^\^Ul^ and a large renormalized Yukawa coupling hk^ 



7-1/27- 



^ip,kn 



1). We 



can therefore neglect in the flow equations all scalar contributions with threshold functions 



involving the large meson masses. This yields the simplified equations [|17| , p!9| for the rescaled 
quantities {d = 4, v^^ = 327r^) 



d_ 
di 

d 



u 



dt 









-4m + (2 + 7]^) ~pu' - ir^l'-o'^i^ph 



Tji^h?' 



N, 



0. 



2 ' 



(58) 



Of course, this approximation is only valid for the initial range of running below fc$ before 
the (dimensionless) renormalized scalar mass squared u'{t,p = 0) approaches zero near the 
chiral symmetry breaking scale. The system (E3) is exactly soluble and we find 



h\t) = Z^\t) 



h] 



u{t,p) 






^^ hjt ' 



N, rt 



Kj 



27r2 



dre-nf^\hi\t)pe 



(59) 



^„2r^ 



Here M/(p) = m(0,p) denotes the effective average potential at the compositeness scale and 
/if is the initial value of h"^ at /c$, i.e. for t = 0. To make the behavior more transparent 
we consider an expansion of the initial value effective potential ui{p) in powers of p around 
p = 



«.(p)^E^p" 



n=0 



n\ 



(60) 



Expanding also /q in eq. (^) in powers of its argument one finds for n > 2 



»^"H^,0) _ 2(n-2)t<M ^ iV, (-l)>-l)! 



/l2«(t) 



hf 



7r2 2"+2(n-2) 



eH(0) 



^2{n-2)t 



(61) 



For decreasing t ^ — 00 the initial values Uj become rapidly unimportant and m("Y^^" 
approaches a fixed point. For n = 2, i.e., for the quartic coupling, one finds 



n(2)(t,0) 


= 1- 


1- 


uf\0) 
hi 


h'{t) 


1- 


Nc l2j. 
8^2 tljt 



(62) 
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McV 


^ 
^ 


fc<J 
McV 




,1/3 

MeV 


rri(fe,j) 
MeV 


m(l GeV) 

MeV 


(VV)(lGeV) 
MeV^ 


Jn 

MeV 


303 


1 


618 


0.0265 


66.8 


14.7 


11.4 


-(186)3 


80.8 


300 





602 


0.026 


66.8 


15.8 


12.0 


-(183)3 


80.2 


310 





585 


0.025 


66.1 


16.9 


12.5 


-(180)3 


80.5 


339 





552 


0.0225 


64.4 


19.5 


13.7 


-(174)3 


81.4 



TABLE I. The table shows the dependence on the constituent quark mass Mq of the input 
parameters k^, m\ /k\ and j as well as some of our "predictions". The phenomenological input 
used here besides Mq is /^r = 92.4 MeV, 171^^ = 135 MeV. The first line corresponds to the values for 
Mq and A/ used in the remainder of this work. The other three lines demonstrate the insensitivity 
of our results with respect to the precise values of these parameters. 



leading to a fixed point value {nP''> /h?):^ = 1. As a consequence of this fixed point behavior 
the system looses all its "memory" on the initial values Uj~ at the compositeness scale 
fc$! Furthermore, the attraction to partial infrared fixed points continues also for the range 
of k where the scalar fluctuations cannot be neglected anymore. However, the initial value 
for the bare dimensionless mass parameter 



<(0) 
h] 



m 



k^ 



A-2 



(63) 



is never negligible. In other words, for hj —>■ oo the IR behavior of the linear quark meson 
model will depend (in addition to the value of the compositeness scale fc$ and the quark mass 
rfi) only on one parameter, ml^. We have numerically verified this feature by starting with 

different values for u] (0). Indeed, the differences in the physical observables were found 
to be small. This IR stability of the flow equations leads to a large degree of predictive 
power! For definiteness we will perform our numerical analysis of the full system of flow 
equations |T^ with the idealized initial value M/(p) = u'j{0)p in the limit hj -^ cxd. It should 
be stressed, though, that deviations from this idealization will lead only to small numerical 
deviations in the IR behavior of the linear quark meson model as long as say hj^l5. 

With this knowledge at hand we may now fix the remaining three parameters of our 



model, fc$, rnl^^ and m by using /^ = 92.4 MeV, the pion mass M^r = 135 MeV and the 
constituent quark mass Mq as phenomenological input. Because of the uncertainty regarding 
the precise value of Mq we give in table | the results for several values of Mg. The first line of 
table I corresponds to the choice of Mq and A/ = m'/(k) which we will use for the forthcoming 
analysis of the model at finite temperature. As argued analytically above the dependence on 
the value of A/ is weak for large enough hj as demonstrated numerically by the second line. 
Moreover, we notice that our results, and in particular the value of j, are rather insensitive 
with respect to the precise value of Mg. It is remarkable that the values for fc$ and rfik^ are 
not very different from those computed in ref. ES 
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IV. HOT QCD AND THE CHIRAL PHASE TRANSITION 
A. Thermal equilibrium and dimensional reduction 

Recall that from the partition function 

Z = Tt {e-^^l = / Dxe-^M (64) 

^ -^ Aperiodic 

all standard thermodjTianiic properties may be determined. The trace over the statistical 
density matrix p = exp(— /5if) with Hamiltonian H and temperature T = l/j3 can be 
expressed as a functional integral |Q, as indicated for a scalar field x with action S by the 
second equation. We supplement Z by a source term J with the substitution S[x] -^ S[x] + 
J Jx in (|6^ . The temperature dependent effective potential U = TT/V for a constant field 
= {x) ) o^^ the Helmholtz free energy, is then 

f/(0;T) = -^lnZj + J0. (65) 

At its minima the effective potential is related to the energy density E/V, the entropy 
density S/V and the pressure p by 

t^lmin = — -T— = -p. (66) 



The term "periodic" in (|64D means that the integration over the field is constrained so 
that 



,s 



$(a;o + (3,x) = $(a;o, x) = ^ <l>,(x)e-^-. , ujf = 2j7tT . (67) 

This is a consequence of the trace operation, which means that the system returns to its 
original state after Euclidean "time" f3. The field can be expanded and we note that the 
lowest Fourier mode for bosonic fields vanishes Uq = 0. For fermionic fields ip the procedure 
is analogous with the important distinction that they are required to be antiperiodic 

7/>(xo + (3,x) = -^(xo, x)=Yl ^,{xy''^'' > < = (2J + 1)^2^ (68) 



as a consequence of the anticommutation property of the Grassmann fields |^ . In contrast 

4 



to the bosonic case the antiperiodicity results in a nonvanishing lowest Fourier mode lo^ 



The extension of flow equations to non-vanishing temperature T is now straightfor- 
ward [^. The (anti-)periodic boundary conditions for (fermionic) bosonic fields in the 



Euclidean time direction leads to the replacement 

in the trace of the flow equation (ESp when represented as a momentum integration, with a 



discrete spectrum for the zero component q'o(j) = ^f ^^i bosons and q'o(j) = ^f for fermions. 
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Hence, for T > a four-dimensional QFT can be interpreted as a three-dimensional model 
with each bosonic or fermionic degree of freedom now coming in an infinite number of copies 
labeled by j G ^ (Matsubara modes). Each mode acquires an additional temperature 
dependent effective mass term q^ij) except for the bosonic zero mode which vanishes identi- 
cally. At high temperature all massive Matsubara modes decouple from the dynamics of the 
system. In this case, one therefore expects to observe an effective three-dimensional theory 
with the bosonic zero mode as the only relevant degree of freedom. One may visualize this 
behavior by noting that for a given characteristic length scale / much larger than the inverse 
temperature jS the compact Euclidean "time" dimension cannot be resolved anymore, as is 
shown in the figure below. This phenomenon is known as "dimensional reduction" . 




_l_ 



P»l 



p«l 



The phenomenon of dimensional reduction can be observed directly from the nonpertur- 
bative flow equations. The replacement ( ]69| ) in ( ^Sj ) manifests itself in the flow equations 
only through a change to T-dependent threshold functions. For instance, the dimensionless 

defined in eq. (p3|) are replaced by 



threshold functions /^(w; r]^ 



T 



n 



5. 



nfl 



Vj'k 



1 U2n- 






d'^-^q ( 1 dRk{q^ 



(2vr; 



d-l 



Z^,k dt J [P{q^) + k^w 



n+l 



(70) 



where q"^ = q^ + q"^ and go = '^'^jT. In the limit k ^ T the sum over Matsubara modes 
approaches the integration over a continuous range of go and we recover the zero temperature 
threshold function /^(i(;;?7$). In the opposite limit fc -C T the massive Matsubara modes 
{I 7^ 0) are suppressed and we expect to find a. d—1 dimensional behavior of l^. In fact, one 
obtains from ([7D| ) 






lUw,T/k;r]^) 



k Vd 



Ut'[ 



for T < A; 



W] rfq,) for T :^ k . 



(71) 



For our choice of the infrared cutoff function Rk, eq. (^q), the temperature dependent 
Matsubara modes in /^(w, T/k; ?7$) are exponentially suppressed for T <^ k. Nevertheless, all 
bosonic threshold functions are proportional to T/k for T ^ k whereas those with fermionic 
contributions vanish in this limit. This behavior is demonstrated in figure ^ where we have 
plotted the quotients l\{w,T /k)/l\{w) and /| {w,T/k)/li (w) of bosonic and fermionic 
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T/k 
FIG. 4. The plot shows the temperature dependence of the bosonic (a) and the fermionic (b) 
threshold functions if {w, T/k) and l\ {w,T/k), respectively, for different values of the dimen- 
sionless mass term w. The solid line corresponds to w = whereas the dotted ones correspond to 
w = 0.1, w = 1 and w = 10 with decreasing size of the dots. For T ^ k the bosonic threshold 
function becomes proportional to T/k whereas the fermionic one tends to zero. In this range the 
theory with properly rescaled variables behaves as a classical three-dimensional theory. 
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threshold functions, respectively. One observes that for k ^ T both threshold functions 
essentially behave as for zero temperature. For growing T or decreasing k this changes as 
more and more Matsubara modes decouple until finally all massive modes are suppressed. 
The bosonic threshold function l\ shows for k <^T the linear dependence on T/k derived in 
eq. (|7TD. In particular, for the bosonic excitations the threshold function for w <^ 1 can be 
approximated with reasonable accuracy by /^(w; //$) for T/k < 0.25 and by {4:T/k)l'}^{w] 77$) 
for T/k > 0.25. The fermionic threshold function l{ tends to zero for k <^T since there 
is no massless fermionic zero mode, i.e. in this limit all fermionic contributions to the flow 
equations are suppressed. On the other hand, the fermions remain quantitatively relevant 
up to T/k ~ 0.6 because of the relatively long tail in figure ^d. The formalism of the average 
action automatically provides the tools for a smooth decoupling of the massive Matsubara 
modes as the momentum scale k is lowered from k ^ T to k <^ T. It therefore allows us 
to directly link the low-T, four-dimensional QFT to the effective three-dimensional high-T 
theory. Whereas for k ^ T the model is most efficiently described in terms of standard four- 
dimensional fields $ a choice of rescaled three-dimensional variables $3 = $/a/T becomes 
better adapted for k <^ T. Accordingly, for high temperatures one will use the rescaled 
dimensionless potential 

k ^ k _ 

U3{t,p3) = -u{t,p) ; p3 = -p . (72) 

For our numerical calculations at non-vanishing temperature we exploit the discussed 
behavior of the threshold functions by using the zero temperature fiow equations in the range 
k > lOT. For smaller values of k we approximate the infinite Matsubara sums (cf. eq. ([70|)) 
by a finite series such that the numerical uncertainty at fc = lOT is better than 10~^. This 
approximation becomes exact in the limit k <^ lOT. 

In section |IIIA| we have considered the relevant fiuctuations that contribute to the fiow 
of Ffc in dependence on the scale k. In thermal equilibrium F^ also depends on the tem- 
perature T and one may ask for the relevance of thermal fiuctuations at a given scale k. 
In particular, for not too high values of T the "initial condition" F^^ for the solution of 
the fiow equations should essentially be independent of temperature. This will allow us to 
fix Ffcj, from phenomenological input at T = and to compute the temperature dependent 
quantities in the infrared (A; -^ 0). We note that the thermal fiuctuations which contribute 
to the r.h.s. of the fiow equation for the meson potential ( |54D are effectively suppressed for 
T^k/A. Clearly for T^/c$/3 temperature effects become important at the compositeness 
scale. We expect the linear quark meson model with a compositeness scale fc$ ~ 600 MeV 
to be a valid description for two fiavor QCD below a temperature of aboutQ 170 MeV. We 
compute the quantities of interest for temperatures T^170MeV by solving numerically the 
T-dependent version of the fiow equations by lowering k from fc$ to zero. For this range of 
temperatures we use the initial values as given in the first line of table |. We observe only a 



^There will be an effective temperature dependence of F^^ induced by the fluctuations of other 
degrees of freedom besides the quarks, the pions and the sigma which are taken into account here. 



We will comment on this issue in section IVD. For realistic three flavor QCD the thermal kaon 



fluctuations will become important for T^lTOMeV. 

32 



minor dependence of our results on the constituent quark mass for the considered range of 
values Mq ~ 300 — 350 MeV. In particular, the value for the critical temperature T^ of the 
model remains almost unaffected by this variation. 



B. High temperature chiral phase transition 

We have pointed out in section | that strong interactions in thermal equilibrium at high 
temperature T differ in important aspects from the well tested vacuum or zero temperature 
properties. A phase transition at some critical temperature T^ or a relatively sharp crossover 



may separate the high and low temperature physics |Q. It was realized early that the 
transition should be closely related to a qualitative change in the chiral condensate according 
to the general observation that spontaneous symmetry breaking tends to be absent in a high 
temperature situation. A series of stimulating contributions PHTII pointed out that for 
sufficiently small up and down quark masses, m^ and md-, and for a sufficiently large mass of 
the strange quark, m^, the chiral transition is expected to belong to the universality class of 
the 0(4) Heisenberg model. It was suggested [0,|r3] ^^at a large correlation length may be 
responsible for important fluctuations or lead to a disoriented chiral condensate. One main 
question we are going to answer using nonperturbative flow equations is: How small ttIu and 
171(1 would have to be in order to see a large correlation length near T^ and if this scenario 
could be realized for realistic values of the current quark masses. 



Figure |^ shows our results |]T^ for the chiral condensate (ipipj as a function of the 
temperature T for various values of the average quark mass m = {niu + md)/2. Curve (a) 
gives the temperature dependence of (ipip) i^i the chiral limit rfi = 0. We flrst consider 
only the lower curve which corresponds to the full result. One observes that the order 
parameter ( V"^) goes continuously (but non-analytically) to zero as T approaches the critical 
temperature in the massless limit Tc = 100.7 MeV. The transition from the phase with 
spontaneous chiral symmetry breaking to the symmetric phase is second order. The curves 
(6), (c) and {d) are for non-vanishing values of the average current quark mass m. The 
transition turns into a smooth crossover. Curve (c) corresponds to mphys or, equivalently, 
m.^{T = 0) = 135 MeV. The transition turns out to be much less dramatic than for m = 0. 
We have also plotted in curve (6) the results for comparably small quark masses ~ 1 MeV, 
i.e. rh = mphys/10, for which the T = value of m^ equals 45 MeV. The crossover is 
considerably sharper but a substantial deviation from the chiral limit remains even for such 
small values of rh. 

For comparison, the upper curves in flgure ^use the universal scaling form of the equation 
of state of the three dimensional 0{4)-sYTainetTic Heisenberg model which will be computed 
explicitly in section [IV C|. We see perfect agreement of both curves in the chiral limit 



for T sufficiently close to T^ which is a manifestation of universality and the phenomenon 
of dimensional reduction. This demonstrates the capability of our method to cover the 
nonanalytic critical behavior and, in particular, to reproduce the critical exponents of the 
0(4)-model (cf. section |I V C|) . Away from the chiral limit we find that the 0(4) universal 
equation of state provides a reasonable approximation for (ipip) in the crossover region 
T = (1.2 - 1.5)Tc. 
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T/T 
FIG. 5. The plot shows the chiral condensate (.Tpi^) as a function of temperature T. Lines (a), 

(b), (c), (d) correspond at zero temperature to m-,^ = 0,45MeV, 135MeV, 230MeV, respectively. 

For each pair of curves the lower one represents the full T-dependence of (iptp) whereas the upper 

one shows for comparison the universal scaling form of the equation of state for the 0(4) Heisenberg 

model. The critical temperature for zero quark mass is Tc = 100.7 MeV. The chiral condensate is 

normalized at a scale /c$ ~ 620 MeV. 



m^/MeV 





45 


135 


230 


Tpc/MeY 


100.7 


-110 


-130 


-150 



TABLE IL The table shows the critical and "pseudocritical" temperatures for various values 
of the zero temperature pion mass. Here Tpc is defined as the inflection point of (tpip) (T). 



In order to facilitate comparison with lattice simulations which are typically performed 
for larger values of m-Tr we also present results for mT^{T = 0) = 230 MeV in curve (d). One 
may define a "pseudocritical temperature" Tpc associated to the smooth crossover as the 
infiection point of (^"ip) (T) as usually done in lattice simulations. Our results for T, 



'-pc 



are 



presented in table ^ for the four different values of m or, equivalently, rriT^iT = 0). The value 
for the pseudocritical temperature for m^ = 230 MeV compares well with the lattice results 
for two flavor QCD (cf. section IVQ) . One should mention, though, that a determination of 



Tpc according to this deflnition is subject to sizeable numerical uncertainties for large pion 
masses as the curve in flgure ^ is almost linear around the inflection point for quite a large 
temperature range. A problematic point in lattice simulations is the extrapolation to realistic 
values of ?ti^ or even to the chiral limit. Our results may serve here as an analytic guide. 
The overall picture shows the approximate validity of the 0(4) scaling behavior over a large 
temperature interval in the vicinity of and above T^ once the (non-universal) amplitudes are 
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T/MeV 
FIG. 6. The plot shows m^ as a function of temperature T for three different values of the average light 

current quark mass m. The solid line corresponds to the realistic value m — mphys whereas the dotted line 

represents the situation without explicit chiral symmetry breaking, i.e., m = 0. The intermediate, dashed 

line assumes m = mphys/ 10. 



properly computed. We point out that the hnk between the universal behavior near Tc and 
zero current quark mass on the one hand and the known physical properties at T = for 
realistic quark masses on the other hand is crucial to obtain all non-universal information 
near T^. 

A second important result is the temperature dependence of the space-like pion cor- 
relation length m~^{T). (We will often call mT^{T) the temperature dependent pion mass 
since it coincides with the physical pion mass for T = 0.) Figure ^ shows mT^{T) and one 
again observes the second order phase transition in the chiral limit m = 0. For T < Tc 
the pions are massless Goldstone bosons whereas for T > Tc they form with the sigma a 
degenerate vector of 0(4) with mass increasing as a function of temperature. For m = 
the behavior for small positive T — Tc is characterized by the critical exponent u, i.e. 
m^{T) = {i+y^Tc{iT -Tc)/TcY and we obtain p = 0.787, ^+ = 0.270. For m > we 
find that m.,^{T) remains almost constant for T^Tc with only a very slight dip for T near 
Tc/2. For T > Tc the correlation length decreases rapidly and for T ^ Tc the precise 
value of m becomes irrelevant. We see that the universal critical behavior near Tc is quite 
smoothly connected to T = 0. The full functional dependence of ?ti^(T, m) allows us to 
compute the overall size of the pion correlation length near the critical temperature and 
we find mT^{Tpc) ~ 1.7m.,^{0) for the realistic value mphys- This correlation length is even 
smaller than the vacuum (T = 0) one and gives no indication for strong fluctuations of 
pions with long wavelength.j] We will discuss the possibility of a tricritical point |]57|,|5|JT^ 



^For a QCD phase transition far from equilibrium long wavelength modes of the pion field can be 
amplified 0,111. 
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FIG. 7. The plot shows the derivative of the meson potential f/(T) with respect to the renor- 

malized field 0/? = (Z<j>p/2)^'^ for different values of T. The first curve on the left corresponds to 

T = 175 MeV. The successive curves to the right differ in temperature by AT = 10 MeV down to 

r = 5MeV. 



with a massless excitation in the two-fiavor case at non-zero baryon number density or for 
three flavors P HTI] , ^ even at vanishing density in section [VE[ We also point out that the 
present investigation for the two flavor case does not take into account a speculative "effec- 
tive restoration" of the axial U a{X) symmetry at high temperature @,^]. We will comment 
on these issues in section [IV V)\ 

In figure |^ we display the derivative of the potential with respect to the renormalized 
field 0R = (Z$p/2)^/^, for different values of T. The curves cover a temperature range 
T = (5 — 175) MeV. The first one to the left corresponds to T = 175 MeV and neighboring 
curves differ in temperature by AT = 10 MeV. One observes only a weak dependence of 
dU{T)/d(j)R on the temperature for T;^60MeV. Evaluated at the minimum of the effective 
potential, 0/j = (Tq, this function connects the renormalized field expectation value with 
mjr(T), the source j and the mesonic wave function renormalization Zi^[T) according to 



dU{T) 



'>R 



o-Qj 



2j 



^R 



.1/2 



Z'^\T) 



AaoiT)ml{T) . 



(73) 



We point out that we have concentrated here only on the meson field dependent part of the 
effective action which is related to chiral symmetry breaking. The meson field independent 
part of the free energy also depends on T and only part of this temperature dependence is 
induced by the scalar and quark fluctuations considered in the present work. Most likely, 
the gluon degrees of freedom cannot be neglected for this purpose. This is the reason why 
we do not give results for "overall quantities" like energy density or pressure as a function 
of T. 

We close this section with a short assessment of the validity of our effective quark meson 
model as an effective description of two flavor QCD at non-vanishing temperature. The 
identification of qualitatively different scale intervals which appear in the context of chiral 
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symmetry breaking, as presented in section [III A| for the zero temperature case, can be gener- 



alized to T 7^ 0: For scales below fc$ there exists a hybrid description in terms of quarks and 
mesons. For k^sB < k^600 MeV chiral symmetry remains unbroken where the symmetry 
breaking scale k^sB(T) decreases with increasing temperature. Also the constituent quark 
mass decreases with T. The running Yukawa coupling depends only mildly on temperature 
for T^120 MeV. fl^ (Only near the critical temperature and for rh = the running is 
extended because of massless pion fluctuations.) On the other hand, for k^AT the effective 
three-dimensional gauge coupling increases faster than at T = leading to an increase of 
AqcniT) with T pDI. As k gets closer to the scale Aqci^{T) it is no longer justified to neglect 



in the quark sector confinement effects which go beyond the dynamics of our present quark 
meson model. Here it is important to note that the quarks remain quantitatively relevant 
for the evolution of the meson degrees of freedom only for scales /c^T/0.6 (cf. figure ^, 
section [IV A]) . In the limit k -C T/0.6 all fermionic Matsubara modes decouple from the 



evolution of the meson potential. Possible sizeable confinement corrections to the meson 
physics may occur if Aqcd(^) becomes larger than the maximum of Mq{T) and T/0.6. This 
is particularly dangerous for small m in a temperature interval around T^. Nevertheless, the 
situation is not dramatically different from the zero temperature case since only a relatively 
small range of k is concerned. We do not expect that the neglected QCD non-localities 
lead to qualitative changes. Quantitative modifications, especially for small rh and |T — T^ 
remain possible. This would only effect the non-universal amplitudes (see sect. [IV C|) . The 



size of these corrections depends on the strength of (non-local) deviations of the quark 
propagator and the Yukawa coupling from the values computed in the quark meson model. 

C. Universal scaling equation of state 



In this section |[T7[ , p!8| we study the linear quark meson model in the vicinity of the 
critical temperature T^ close to the chiral limit m = 0. In this region we find that the sigma 
mass m~^ is much larger than the inverse temperature T~^, and one observes an effectively 
three-dimensional behavior of the high temperature quantum field theory. We also note 
that the fermions are no longer present in the dimensionally reduced system as has been 



discussed in section P!VA[ We therefore have to deal with a purely bosonic 0(4)-symmetric 
linear sigma model. At the phase transition the correlation length becomes infinite and 
the effective three-dimensional theory is dominated by classical statistical fluctuations. In 
particular, the critical exponents which describe the singular behavior of various quantities 
near the second order phase transition are those of the corresponding classical system. 

Many properties of this system are universal, i.e. they only depend on its symmetry 
(0(4)), the dimensionality of space (three) and its degrees of freedom (four real scalar com- 
ponents). Universality means that the long-range properties of the system do not depend on 
the details of the specific model like its short distance interactions. Nevertheless, important 
properties as the value of the critical temperature are non-universal. We emphasize that 
although we have to deal with an effectively three-dimensional bosonic theory, the non- 
universal properties of the system crucially depend on the details of the four-dimensional 
theory and, in particular, on the fermions. 

Our aim is a computation of the scaling form of the equation of state which relates for 
arbitrary T near T^. the derivative of the free energy or effective potential U to the average 
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current quark mass rh. At the critical temperature and in the chiral hmit there is no scale 
present in the theory. In the vicinity of T^ and for small enough rh one therefore expects 
a scaling behavior of the dimensionless average potential Uk = k~'^Uk as a function of the 
rescaled field p = Z^^k^^'^p |]6T| . |6^ . (See also the review in ||63[| .) 



There are only two independent scales close to the transition point which can be related 
to the deviation from the critical temperature, T — Tc, and to the explicit symmetry breaking 
by a nonvanishing quark mass rh. As a consequence, the properly rescaled potential can 
only depend on one scaling variable. A possible choice is the Widom scaling variable 

(T-T^/n 
{2ao/T; 



-^^:^4^- (74) 



Here f3 is the critical exponent of the order parameter Wq in the chiral limit rh = (see 
equation ([78|)). With U'{p = 2ctq) = j/(2ao) the Widom scaling form of the equation of 
state reads [Q 



' ^^"'^'m (75) 



where the exponent 6 is related to the behavior of the order parameter according to (pOl). The 
equation of state (|75D is written for convenience directly in terms of four-dimensional quan- 
tities. They are related to the corresponding effective variables of the three-dimensional 
theory by appropriate powers of Tc. The source j is determined by the average current 
quark mass m as j = 2m\^^rh. The mass term at the compositeness scale, rn\^, also re- 
lates the chiral condensate to the order parameter according to yl)il>j = — 2m|_^(ao — rh). 
The critical temperature of the linear quark meson model was found in section |1V lj| to be 
Tc = 100.7 MeV. 

The scaling function / is universal up to the model specific normalization of x and itself. 
Accordingly, all models in the same universality class can be related by a rescaling of Wq 
and T — Tc. The non-universal normalizations for the quark meson model discussed here 
are defined according to 

f{Q)=D , /(-5-V/3)=o. (76) 

We find D = 1.82 ■ 10~^, B = 7.41 and our result for (3 is given in table |T|. Apart from the 
immediate vicinity of the zero of / (x) we find the following two parameter fit for the scaling 
function, 

/fit(x) = 1.816 ■10^^(1 + 136.1 x)2(l + 160.9 ^x)"^ 

(77) 



1 + 160.9(0.9446^^) 



-l/(7-2-A) 2;)7-2-A 



to reproduce the numerical results for / and df /dx at the 1 — 2% level with = 0.625 (0.656), 
A = —0.490 (—0.550) for x > (x < 0) and 7 as given in table 0. The universal properties 
of the scaling function can be compared with results obtained by other methods for the 
three-dimensional 0(4) Heisenberg model. In figure ^ we display our results along with 



those obtained from lattice Monte Carlo simulation |65], second order epsilon expansion 



and mean field theory. We observe a good agreement of average action, lattice and 
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FIG. 8. The figure shows a comparison of our results'^ denoted by "average action", with results 
of other methods for the scaling function of the three-dimensional 0(4) Heisenberg model. We 
have labeled the axes for convenience in terms of the expectation value ctq and the source j of 
the corresponding four-dimensional theory. The constants B and D specify the non-universal 
amplitudes of the model (cf . eq. ^) . The curve labeled by "MC" represents a fit to lattice Monte 
Carlo data. The second order epsilon expansion [BH] and mean field results are denoted by "e" and 



"mf", respectively. Apart from our results the curves are taken from ref. |65]. 
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epsilon expansion results within a few per cent for T < Tc. Above Tc the average action and 
the lattice curve go quite close to each other with a substantial deviation from the epsilon 
expansion and mean field scaling function. (We note that the question of a better agreement 
of the curves for T < T^ or T > T^ depends on the chosen non-universal normalization 
conditions for x and /, eq. ([76|).) 

Before we use the scaling function f{x) to discuss the general temperature and quark 
mass dependent case, we consider the limits T = Tc and m = 0, respectively. In these 
limits the behavior of the various quantities is determined solely by critical amplitudes and 
exponents. In the spontaneously broken phase (T < T^) and in the chiral limit we observe 
that the renormalized and unrenormalized order parameters scale according to 



rji \ / \ rji 

; ^ (78) 

2ao(T) fTc-T\^ 

Jd 



Tc V Tc 



respectively, with E = 0.814 and the value of B given above. In the symmetric phase the 

V 



renormalized mass m = m.,^ = rria and the unrenormalized mass m = Zj m behave as 



miT) ,^^-ifT-T, 



(79) 
where ^^ = 0.270, C^ = 2.79. For T = Tc and non-vanishing current quark mass we have 

with the value of D given above. 

Though the five amplitudes E, B, ^+, C^ and D are not universal there are ratios of 
amplitudes which are invariant under a rescaling of a^ and T — Tc. Our results for the 
universal amplitude ratios are 

R^ = C+DB^-^ = 1.02 , 

R^ = (^+)/5/-£)i/(5+i)5 = 0.852 , (81) 

^+E = 0.220 . 

Those for the critical exponents are given in table fT|. For comparison table |T| also gives 
the results from perturbation series at fixed dimension to seven-loop order |]67| , |69| as well 
as lattice Monte Carlo results |^ which have been used for the lattice form of the scaling 



function in figure En There are only two independent amplitudes and critical exponents, 



^See also ref. [^ and references therein for a calculation of critical exponents using similar 
methods as in this work. 
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V 


7 


b 


/? 


7? 


average action 


0.787 


1.548 


4.80 


0.407 


0.0344 


FD 


0.73(2) 


1.44(4) 


4.82(5) 


0.38(1) 


0.03(1) 


MC 


0.7479(90) 


1.477(18) 


4.851(22) 


0.3836(46) 


0.0254(38) 



TABLE III. The table shows the critical exponents corresponding to the three-dimensional 
0(4)-Heisenberg model. Our results are denoted by "average action" whereas "FD" labels the 



exponents obtained from perturbation series at fixed dimension to seven loops |67]. The bottom 



line contains lattice Monte Carlo results |68|. 



respectively. They are related by the usual scaling relations of the three-dimensional scalar 
0(A^)-model |^9| which we have explicitly verified by the independent calculation of our 
exponents. 

We turn to the discussion of the scaling behavior of the chiral condensate v^w) for the 
general case of a temperature and quark mass dependence. In figure ^ we have displayed 
our results for the scaling equation of state in terms of the chiral condensate 



'\\)'\\) 



—my 



-Tr 






i/s 



+ J 



(82) 



as a function of T/T^. = 1 + x{j/T'^f{x)Y^^^ for different quark masses or, equivalently, 
different values of j. The curves shown in figure ^ correspond to quark masses m = 0, 
^ = ^phys/10, rh = rhphys and rh = 3.5mphys or, equivalently, to zero temperature pion 
masses m.,^ = 0, m^r = 45 MeV, 771.,^ = 135 MeV and m^ = 230 MeV, respectively. The scaling 
form (^) for the chiral condensate is exact only in the limit T ^ T^, j ^ 0. It is interesting 



to find the range of temperatures and quark masses for which (ipip) approximately shows the 
scaling behavior (^2]). This can be infered from a comparison (see figure ^ with our full non- 
universal solution for the T and j dependence of ("ip^) as described in section [IV B|. For m^ = 



one observes approximate scaling behavior for temperatures T ^ 90 MeV. This situation 
persists up to a pion mass of rrij^ = 45 MeV. Even for the realistic case, m^ = 135 MeV, 
and to a somewhat lesser extent for ?Ti,r = 230 MeV the scaling curve reasonably reflects the 
physical behavior for T^T^. For temperatures below Tc-, however, the zero temperature mass 
scales become important and the scaling arguments leading to universality break down. 

The above comparison may help to shed some light on the use of universality arguments 
away from the critical temperature and the chiral limit. One observes that for temperatures 
above Tc the scaling assumption leads to quantitatively reasonable results even for a pion 
mass almost twice as large as the physical value. This in turn has been used for two 
flavor lattice QCD as theoretical input to guide extrapolation of results to light current 
quark masses. From simulations based on a range of pion masses 0.3^m7r/mp^0.7 and 
temperatures < T^250MeV a "pseudocritical temperature" of approximately 140 MeV 
with a weak quark mass dependence is reported |]7T|. Here the "pseudocritical temperature" 



Tpc is defined as the inflection point of (V'V') as a function of temperature. For comparison 



41 



with lattice data we have displayed in figure ^ the temperature dependence of the chiral 
condensate for a pion mass 771^^ = 230 MeV. From the free energy of the linear quark meson 
model we obtain in this case a pseudocritical temperature of about 150 MeV in reasonable 
agreement with lattice results ||7T1 , |72| . In contrast, for the critical temperature in the chiral 
limit we obtain Tc = 100.7 MeV. This value is considerably smaller than the lattice results 



of about (140 — 150) MeV obtained by extrapolating to zero quark mass in ref. [^. We 
point out that for pion masses as large as 230 MeV the condensate (i/ji/j) (T) is almost linear 
around the inflection point for quite a large range of temperature. This makes a precise 
determination of T^, somewhat difficult. Furthermore, figure |^ shows that the scaling form 
of (i/ji/j) (T) underestimates the slope of the physical curve. Used as a fit with Tc as a 
parameter this can lead to an overestimate of the pseudocritical temperature in the chiral 
limit. 

The linear quark meson model exhibits a second order phase transition for two quark 
fiavors in the chiral limit. As a consequence the model predicts a scaling behavior near the 
critical temperature and the chiral limit which can, in principle, be tested in lattice simu- 
lations. For the quark masses used in the present lattice studies the order and universality 
class of the transition in two flavor QCD remain a partially open question. Though there 
are results from the lattice giving support for critical scaling there are also simulations with 



two flavors that reveal signiflcant flnite size effects and problems with 0(4) scaling [^. 



D. Additional degrees of freedom 

So far we have investigated the chiral phase transition of QCD as described by the linear 
0(4)-model containing the three pions and the sigma resonance as well as the up and down 
quarks as degrees of freedom. Of course, it is clear that the spectrum of QCD is much richer 
than the states incorporated in our model. It is therefore important to ask to what extent 
the neglected degrees of freedom like the strange quark, strange (pseudo)scalar mesons, 
(axial) vector mesons, baryons, etc., might be important for the chiral dynamics of QCD. 
Before doing so it is instructive to flrst look into the opposite direction and investigate the 
difference between the linear quark meson model described here and chiral perturbation 



theory based on the non-linear sigma model [22|. In some sense, chiral perturbation theory 



is the minimal model of chiral symmetry breaking containing only the Goldstone degrees 
of freedom]^. By construction it is therefore only valid in the spontaneously broken phase. 
For small temperatures (and momentum scales) chiral perturbation theory is expected to 
give a reliable description whereas for high temperatures or close to Tc one expects sizeable 
corrections. 

From |7^ we infer the three-loop result for the temperature dependence of the chiral 
condensate in the chiral limit for A^ light flavors 



^ For vanishing temperature it has been demonstrated [ 73 , 7J] that the results of chiral perturbation 
theory can be reproduced within the linear meson model once certain higher dimensional operators 
in its effective action are taken into account for the three flavor case. 
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FIG. 9. The plot displays the chh'al condensate (ipip) as a function of T/fir ' . The solid line 
corresponds to our results for vanishing average current quark mass m = whereas the dashed 
line shows the corresponding three-loop chiral perturbation theory result for Fi = 470 MeV. 
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The scale Fi can be determined from the D-wave isospin zero tttt scattering length and is 
given by Fi = (470 ± 100) MeV. The constant Fq is (in the chiral limit) identical to the pion 
decay constant Fq = /^°^ = 80.8 MeV (cf. table |). In figure ^ we have plotted the chiral 



condensate as a function of T/Fq for both, chiral perturbation theory according to (^) and 
for the linear quark meson model. As expected the agreement for small T is very good. 
Nevertheless, small numerical deviations are manifest even for T -C Tc due to quark and 
sigma meson loop contributions. We observe for larger values of T, say for T ^ 0.8/^°-*, that 
the deviations become significant. Chiral perturbation theory is not expected to correctly 
reproduce the critical behavior of the system near its second order phase transition. 

Within the language of chiral perturbation theory the neglected effects of thermal quark 
fluctuations may be described by an effective temperature dependence of the parameter 
Fq{T). We notice that the temperature at which these corrections become important equals 
approximately one third of the constituent quark mass Mq{T) or the sigma mass mfj{T), 
respectively, in perfect agreement with figure §. As suggested by this figure the onset of 
the effects from thermal fluctuations of heavy particles with a T-dependent mass mniT) 
is rather sudden for T^mHiT)/3. These considerations also apply to our two flavor quark 
meson model. Within full QCD we expect temperature dependent initial values at k<^. 

The dominant contribution to the temperature dependence of the initial values presum- 
ably arises from the influence of the mesons containing strange quarks as well as the strange 
quark itself. Here the quantity ml^ seems to be the most important one. (The temperature 
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dependence of higher couphngs hke A(T) is not very relevant if the IR attractive behavior 
remains vahd, i.e. if 2^$,^^ remains small for the range of temperatures considered. We ne- 
glect a possible T-dependence of the current quark mass m.) In particular, for three flavors 
the potential Uk^ contains a term 

- 2'^k^ (^^t '^ + ^^t '^^) = -'i'k^Vs'^uu'^dd + ■■■ (84) 

which reflects the axial Ua{^) anomaly. It yields a contribution to the effective mass term 
proportional to the expectation value ($ss) = V^s, i-e. 

^^L = -2'^k^Vs ■ (85) 

Both, I7fc_j^ and (ps, depend on T. We expect these corrections to become relevant only for 
temperatures exceeding m;^(T)/3 or Ms{T)/3. We note that the temperature dependent 
kaon and strange quark masses, mK{T) and Ms{T), respectively, may be somewhat different 
from their zero temperature values but we do not expect them to be much smaller. A 
typical value for these scales is around 500 MeV. Correspondingly, the thermal fluctuations 
neglected in our model should become important for T^170 MeV. It is even conceivable 
that a discontinuity appears in (ps{T) for sufficiently high T (say T ~ 170 MeV). This would 
be reflected by a discontinuity in the initial values of the 0(4)-model leading to a first 
order transition within this model. Obviously, these questions should be addressed in the 
framework of the three fiavor ^^/^(S) x SUii{3) quark meson model. Work in this direction 
is in progress. 

We note that the temperature dependence oiT'(T)(ps{T) is closely related to the question 

of an effective high temperature restoration of the axial f/A(l) symmetry PJ59[|. The r]' 

mass term is directly proportional to this combination, m'i,{T) — m^{T) ~ ^V{T)ips{T) 

76| . Approximate Ua{X) restoration would occur if (ps(T) or T'iT) would decrease sizeable 



for large T. For realistic QCD this question should be addressed by a three flavor study. 
Within two flavor QCD the combination V^ips is replaced by an effective anomalous mass 
term 17^ . We add that this question has also been studied within full two flavor QCD in 



lattice simulations [^ but no flnal conclusion can be drawn yet. 

To summarize, we have found that the effective two flavor quark meson model presumably 
gives a good description of the temperature effects in two flavor QCD for a temperature range 
T^170MeV. Its reliability should be best for low temperature where our results agree 
with chiral perturbation theory. However, the range of validity is considerably extended as 
compared to chiral perturbation theory and includes, in particular, the critical behavior of 
the second order phase transition in the chiral limit. The method using nonperturbative 
flow equations within the linear quark meson model may also help to shed some light on the 
remaining pressing questions at high temperature, like the nature of the phase transition for 
realistic values of the strange quark mass. 
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V. HIGH BARYON NUMBER DENSITY 

Over the past years, considerable progress has been achieved in our understanding of 
high temperature QCD, where simulations on the lattice and universality arguments played 
an essential role. The results of the renormalization group approach to the effective quark 
meson model provides the link between the low temperature chiral perturbation theory 



domain of validity and the high temperature domain of critical phenomena ^j. On the 
other hand, our knowledge of the high density properties of strongly interacting matter is 
rudimentary so far. There are severe problems to use standard simulation algorithms at 



nonzero chemical potential on the lattice because of a complex fermion determinant |7^. In 
this section we will try to get some insight into matter at high density using nonperaturbative 
flow equations 



A. Cold dense matter 

We will consider here a few general aspects of QCD at nonzero baryon number density 
and vanishing temperature. For a more detailed recent review see ref. |12|. It is instructive 
to consider for a moment a free theory of fermions with mass m carrying one unit of baryon 
charge where the associated chemical potential is /x. When yU > the ground state is 
the Fermi sphere with radius pp = \fj^^^^rr? and therefore n(/i) = (//^ — m^)^/^/(37r^). 
For n < m the density vanishes identically and n represents an order parameter which 
distinguishes the two phases. One clearly observes the nonanalytic behavior at /x = m 
which denotes the "onset" value for nonzero density. 

Nuclear liquid gas transition. — What is the "onset" value for nonzero density in QCD? 
The low density properties of QCD may be inferred from what is known empirically about 
bulk nuclear matter. Here nuclear matter denotes a uniform, isospin symmetric large sample 
of matter. Uniform infinite symmetric matter is an idealization whose properties, however, 
are related to those of finite nuclei in the liquid droplet model of the nucleus. Based on 
information about finite nuclei one finds an energy per baryon of m]\f — 16 MeV where 
mjsf ~ rUn ~ nip is the nucleon mass and 16 MeV is the binding energy per nucleon. In 
QCD one expects the density to jump at /Xnuc — f^N — 16 MeV from zero to saturation 
density uq at which the pressure is zero. A schematic plot is shown in figure |l^. For small 
temperature this corresponds to the transition between a gas of nucleons and nuclear matter. 
The situation changes once other interactions, in particular electromagnetic interactions, are 
taken into account. Including electromagnetism and adding electrons for charge neutrality 
the "onset" density is that of the lowest energy state of hadronic matter - ordinary iron 
with a density about 10~^^ times nuclear matter density. 

Chiral phase transition — We have reasonable expectations, but much less evidence, for 
the behavior of QCD at densities much larger than nuclear matter density. One may suppose 
as a starting point that at very high density quarks behave nearly freely and form large Fermi 
surfaces. The leading interactions for particle-hole excitations near the Fermi surface are 
elastic scatterings among these particles and holes, since other possibihties are blocked by 
the Pauli exclusion principle. Because the momenta are large, scattering processes typically 
involve large momentum transfers and are therefore indeed characterized by small couplings 
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M-nuc M-qm 

FIG. 10. Schematic plot of the baryon number density in QCD as a function of the chemical 

potential for zero temperature. 



due to asymptotic freedom. (This self-consistent argument can be worked out to include 
also the neglected "small angle" scatterings.) 

What can we expect at a chemical potential /i ~ fi^m (cf. figure 10) associated with the 
quark hadron transition, or at densities of about two to ten times nuclear matter density? 
In vacuum chiral symmetry breaking is caused by a condensate of quark-ant iquark pairs 
with zero net momentum. Since at nonzero density all particle states up to the Fermi 
momentum pp are occupied, there is no possibility for correlated pairs at low momenta, and 
at high momentum the quark-antiquark pairs have large energy (cf. also figure |l^ in section 
|V b\ ) . At sufficiently high density we therefore expect the chiral condensate {ipi') to drop 
significantly. On the other hand there are attractive quark-quark channels. A particularly 
interesting possibility is the formation of Cooper pairs of quarks at high density which 
condense. We will discuss this striking feature in section [VF| . First we will concentrate on 
the fundamental question for QCD about the nature of the transition that restores chiral 
symmetry. 



B. The quark meson model at nonzero density 

We extend our discussion of the linear quark meson model to nonzero baryon number 
density ||2^. The model can be viewed as a generalization of the Nambu-Jona-Lasinio 
model for QCD ||5ll. Calculations at nonzero baryon density are typically based on a mean 

that the order of the chiral phase transition 



field approximation and it was claimed ||51| J 
is ambiguous. Recently, the order of this high density transition has raised considerable 
interest. One can argue on general grounds P,p!2[] that if two-flavor QCD has a second order 
transition (crossover) at high temperatures and a first order transition at high densities, 
then there exists a tricritical point (critical endpoint) with long-range correlations in the 
phase diagram. As pointed out in section |, the physics around this point is governed by 
universality and may allow for distinctive signatures in heavy ion collisions [^,|T2|,|T3|. We 
will discuss this important issue in section [VE| . 

We will observe that in a proper treatment beyond mean field theory the order of the 
chiral phase transition can be fixed within the models under investigation. For this purpose, 
a crucial observation is the strong attraction of the flow to partial infrared fixed points 
discussed in section p^llB] |T9|JT7[] . The two remaining relevant or marginal parameters can 
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be fixed by the phenomenological values of /^ and tlie constituent quark mass. For two 
massless quark flavors we find that chiral symmetry restoration occurs via a first order 
transition between a phase with low baryon density and a high density phase. We emphasize, 
nevertheless, that the linear quark meson model captures the low density properties of QCD 
only incompletely since the effects of confinement are not included. In particular, for a 
discussion of the liquid-gas nuclear transition the inclusion of nucleon degrees of freedom 
seems mandatory |7^. On the other hand this model is expected to provide a reasonable 
description of the high density properties of QCD. 

In quantum field theory the effects of a non-vanishing baryon density in thermal equi- 
librium or the vacuum are described by adding to the classical action a term proportional 
to the chemical potential n, 

A^S = tfiJ2 fj^ dxo J d^x ^^^% = -3^B . (86) 

For quarks the sum is over A'"c colors and Np flavors. With our conventions /x corresponds 
to the chemical potential for the quark number density. The baryon number density n can 
be obtained from the /x-dependence of the Euclidean effective action F, evaluated at its 
minimum for fixed temperature T and volume[3 V 

{B) 1 d F^i,T 



n = 



(87) 

T,V 



V 3dfi V 

We will restrict our discussion here to two massless quarks. A more realistic treatment 
would have to include the small up and down quark masses and the finite, but much heavier, 
strange quark mass. Though their inclusion will change certain quantitative estimates, they 
most likely do not change the qualitative outcome of the investigation as is discussed below. 
The linear quark meson model is defined in section ^. At nonzero temperature T and 
chemical potential /i our ansatz for F^ reads 



Tt = £'^ dx" j cPxUi'(',^^^, + n'')i'a + hti,' 



ii^*.'-i^(**).' 



^b 



+Z$,fc9^$:,a'^$"'' + Uk{-p\ /x, T) I . (88) 

A nonzero chemical potential yU to lowest order results in the term ~ i^ijj •y'^ipa appearing 
on the right hand side of (^). Our approximation neglects the dependence of Z^^k and 
hk on /i and T. We also neglect a possible difference in normalization of the quark kinetic 
term and the baryon number current. The form of the effective action at the compositeness 
scale, Ffc^J-?/;, $], serves as an initial value for the renormalization group flow of Tk[ip, $] for 
k < kq,. We will consider here the case that Z^ ^^ ^ 1. The limiting case Z<^^k^ = can be 
considered as a solution of the corresponding Nambu-Jona-Lasinio model where the effective 
four-fermion interaction has been eliminated by the introduction of auxiliary meson fields 
(cf. section pi]) . 



-^^More precisely, B counts the number of baryons minus antibaryons. For T — > the factor T/V 
is simply the inverse volume of four dimensional Euclidean space. 
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C. Renormalization group flow 

1. Flow equation for the effective potential 

The dependence on the infrared cutoff scale k of the effective action r^ is given by the 
exact flow equation (^). We employ the same infrared cutoff function for the bosonic fields 
RkB (il) as in the previous sections. A few comments about the fermionic cutoff function 
are in place. The infrared cutoff function for the fermions RkF should be consistent with 
chiral symmetries. This can be achieved if RkF has the same Lorentz structure as the kinetic 
term for free fermions [^ . In presence of a chemical potential /i we use 

RkF = (Yq^ + it^l^kF . (89) 

The effective squared inverse fermionic propagator is then of the form 

PkF = [(go + if^f + g'](l + rkF? 

= (go + ^/i)' + f + A;'e(4 - (go + ^/i)' - f) , (90) 

where the second line defines VkF and one observes that the fermionic infrared cutoff acts as 
an additional mass-like term ~ k^^ Here the B-function can be thought of as the limit of 
some suitably regularized function, e.g. 9^ = [exp{(go + z/i)^ + g ^ — fc|}/e + 1]~^- 

We compute the flow equation for the effective potential Uk from equation (|i3| ) using the 
ansatz ( ^8]) for r^ and we introduce a renormalized field p = Zi^^kp and Yukawa coupling 
hk = Z^ i! hk- In complete analogy to section |III B| we find that the flow equation for f/j 



k 



obtains contributions from bosonic and fermionic fluctuations, respectively. 



-f/.Mp;T,/.) = lTE/ 



(27r)3Zc,,fe dk ]PkB{q^) + U[{p;T,^i) 



^ PkB (g2) + U'kip; T, p) + 2pU';{p- T,p)]' ^^^^ 

9 TT ( T \ &Mi-\^ 7 ^^g ke{kl-{qo + ipf-q^) 

— UkF{p;T,p) = -8NcT}^ / , ,3 ,, , . ^2 , ^^ , , 2 /o • ^2) 

ok ^J [2Txy'' PkF[{qQ + ipY + q^) + Kp2 



^^The exponential form (^) of the cutoff function R^q renders the first term on the right hand 
side of ( |45| ) both infrared and ultraviolet finite. No need for an additional ultraviolet regularization 
arises in this case. This is replaced by the necessary specification of an initial value T^^ at the 
scale fc$. Here /c$ is associated with a physical ultraviolet cutoff in the sense that effectively 
all fluctuations with q^ > k\ are already included in T^^. A similar property for the fermionic 
contribution is achieved by the G-function in (pO|). We note that in the previous sections an 
exponential form of the fermionic infrared cutoff function R^f was used. At nonzero density the 
mass-like IR cutoff simplifies the computations considerably because of the trivial momentum 
dependence. 
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Here q"^ = Qq + (p with q^in) = 2n7cT for bosons, qo{n) = (2n + 1)ttT for fermions {n G Z) 
and Nc = 3 denotes the number of colors. The scale dependent propagators on the right 
hand side contain the momentum dependent pieces PkB = q"^ + Z^\Rk{q'^) and PkF given 
by (^) as well as mass terms. The only explicit dependence on the chemical potential n 
appears in the fermionic contribution (^) to the flow equation for Uk- It is instructive to 
perform the summation of the Matsubara modes explicitly for the fermionic part. Since 
the flow equations only involve one momentum integration, this can be easily done with 
standard techniques using contour integrals. (See e.g. IQ chapter 3.) One finds 



'-U,Ap;T,,)^-8N. I ^ '^^frS\. +^N. ^ ''' 



dk '^^^^' '"^ \l {2nyq^ + k^ + hlp/2' ^^ (2vr)3 ^^ + p + ;,2^/2 

(93) 



X 



exp 



q^ + k^ + hlp/2-fi)/T +1 exp {Jq^ + k^ + hlp/2 + fi)/T 



where, for simplicity, we sent A;$ — > oo in the /i, T dependent second integral. This is justified 
by the fact that in the /z, T dependent part the high momentum modes are exponentially 
suppressed. 

For comparison, we note that within the present approach one obtains standard mean 
field theory results for the free energy if the meson fluctuations are neglected, dUks/dk = 0, 



and the Yukawa coupling is kept constant, hk = h m. (|93|). The remaining flow equation 
for the fermionic contribution could then easily be integrated with the (mean field) initial 
condition Uk^{p) = rn\^p. In the following we will concentrate on the case of vanishing tem- 
perature. We find (see below) that a mean field treatment yields relatively good estimates 
only for the /^-dependent part of the free energy U{p] /x) — t/(p; 0). On the other hand, mean 
field theory does not give a very reliable description of the vacuum properties which are 
important for a determination of the order of the phase transition at /x 7^ 0. 

2. Zero temperature physics 

In the limit of vanishing temperature one expects and observes a non-analytic behavior 
of the /i-dependent integrand of the fermionic contribution (^) to the flow equation for U^ 
because of the formation of Fermi surfaces. Indeed, the explicit //-dependence of the flow 
equation reduces to a step function 

9rr ...,.^ o.r 7 A ke{kl - q') 



-UkFip; f^) = - 8N, J 



dk '^^ ''^' ^' ' J (27r)4 g2 + A;2 + hlp/2 

— 00 

+ 4iV, / ,f ^3 , ^ = e(p- Jq^ + k^ + hlp/2) . (94) 

The quark chemical potential p enters the bosonic part (^TJ) of the flow equation only 
imphcitly through the meson mass terms Uj^^p; p) and U^^p] p) + 2pf/^'(p;/x) for the pions 
and the cr-meson, respectively. For scales k > p the B-function in (0) vanishes identically 
and there is no distinction between the vacuum evolution and the p ^ evolution. This 
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is due to the fact that our infrared cutoff adds to the effective quark mass (fc^ + /i|p/2)^/^. 
For a chemical potential smaller than this effective mass the "density" —dUk/dn vanishes 
whereas for larger /x one can view /i = [^^(/i, k, p) + k^ + /i|p/2]^/^ as an effective Fermi 
energy for given k and p. A small infrared cutoff k removes the fluctuations with momenta 
in a shell close to the physical Fermi surfaceQ p'^ — h\p/2 — k"^ < q^ < p'^ — h\^Qp/2. Our 
flow equation realizes the general idea [Q that for p y^ the lowering of the infrared cutoff 
/c — > should correspond to an approach to the physical Fermi surface. For a computation 
of the meson effective potential the approach to the Fermi surface in (|9^) proceeds from 
below and for large k the effects of the Fermi surface are absent. By lowering k one "fills 
the Fermi sea". 

As discussed in section [111B| the observed fixed point behavior in the symmetric regime 
allows us to fix the model by only two phenomenological input parameters and we use /^r = 
92.4 MeV and 300MeV^Mq^350MeV. The results for the evolution in vacuum described 



in section |ITT| [|T^,|T^ show that for scales not much smaller than k^ ~ 600 MeV chiral 
symmetry remains unbroken. This holds down to a scale of about k^sB — 400 MeV at which 
the meson potential Uk{p) develops a minimum at po.fc > thus breaking chiral symmetry 
spontaneously. Below the chiral symmetry breaking scale running couplings are no longer 
governed by the partial fixed point. In particular, for k^k^sB the Yukawa coupling h^ and 
the meson wave function renormalization Z$ ^ depend only weakly on k and approach their 
infrared values. At p 7^ we will follow the evolution from k = k^sB to k = and neglect 
the fc-dependence of hk and Z$^fc in this range. According to the above discussion the initial 
value Uk^gg is /i-independent for p < k^sB- We solve the flow equation ( p6D with (pTD, 
(0) numerically as a partial differential equation for the potential depending on the two 
variables p and k for given p. Nonzero current quark masses, which are neglected in our 
approximation, result in a pion mass threshold and would effectively stop the renormalization 
group flow of renormalized couplings at a scale around m^. To mimic this effect one may 
stop the evolution by hand at kj ~ m^. We observe that our results are very insensitive 
to such a procedure which can be understood from the fact that a small infrared cutoff 
- induced by kf or by the nonzero current quark masses - plays only a minor role for a 
sufficiently strong first order transition.[^ 

In the fermionic part (P^) of the flow equation the vacuum and the p-dependent term 



contribute with opposite signs. This cancelation of quark fluctuations with momenta below 

^^If one neglects the mesonic fluctuations one can perform the ^-integration of the flow equa- 
tion (p^ ) in the limit of a /^-independent Yukawa coupling. One recovers (for /c|, ^ fc^ + h'^p/2, p'^) 
mean field theory results except for a shift in the mass, h?p/2 -^ h? p/2 + k"^ , and the fact that 
modes within a shell of three-momenta p^ — h^p/2 — k"^ < q"^ < p^ — h^p/2 are not yet included. 
Because of the mass shift the cutoff k also suppresses the modes with q^ < k"^. 

^^For the results presented in the next section we use kf = 100 MeV. The flow of the potential U^ 
around its minima and for its outer convex part stabilizes already for k somewhat larger than kf. 
Fluctuations on larger length scales lead to a flattening of the barrier between the minima. The 
approach to convexity is not relevant for the present discussion. 
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the Fermi surface is crucial for the restoration of chiral symmetry at high density]^. In 
vacuum, spontaneous chiral symmetry breaking is induced in our model by quark fluctuations 
which drive the scalar mass term Ul^{p = 0) from positive to negative values at the scale 
k = k^sB- (Meson fluctuations have the tendency to restore chiral symmetry because of 
the opposite relative sign, cf. (^Sf).) As the chemical potential becomes larger than the 
effective mass (fc^ + h\p/2y/'^ quark fluctuations with momenta smaller than qp{n,k,p) = 
fi^ — k"^ — h\p/2 are suppressed. Since qp is monotonically decreasing with p for given p and 
k the origin of the effective potential is particularly affected. We will see in the next section 
that for large enough p this leads to a second minimum of Uk=o{p; p) ai p = and a chiral 
symmetry restoring flrst order transition. 

D. High density chiral phase transition 



In vacuum or at zero density the effective potential U{a), a = J p/2, has its minimum 
at a nonvanishing value ctq = f-w/'^ corresponding to spontaneously broken chiral symmetry. 
As the quark chemical potential p increases, U can develop different local minima. The 
lowest minimum corresponds to the state of lowest free energy and is favored. In flgure ^ 
we plot the free energy as a function of a for different values of the chemical potential 
p = 322.6, 324.0, 325.2 MeV. Here the effective constituent quark mass is M, = 316.2 MeV. 
We observe that for p < Mq the potential at its minimum does not change with p. The 
quark number density is 

'-'r' |min 

and we conclude that the corresponding phase has zero density. In contrast, for a chemical 
potential larger than Mg we flnd a low density phase where chiral symmetry is still broken. 



The quark number density as a function of p is shown in flgure |T2[ One clearly observes 
the non-analytic behavior at p = Mg which denotes the "onset" value for nonzero density. 
From flgure |TI| one also notices the appearance of an additional local minimum at the origin 
of U. As the pressure p = —U increases in the low density phase with increasing /i, a 
critical value pc is reached at which there are two degenerate potential minima. Before p 
can increase any further the system undergoes a flrst order phase transition at which two 
phases have equal pressure and can coexist. In the high density phase chiral symmetry 
is restored as can be seen from the vanishing order parameter for p > p^. We note that 
the relevant scale for the flrst order transition is Mg. For this reason we have scaled our 
results for dimensionful quantities in units of Mg. For the class of quark meson models 
considered here (with Mg/ f^^ in a realistic range around 3-4) the flrst order nature of the 
high density transition has been clearly established. In particular, these models comprise 
the corresponding Nambu-Jona-Lasinio models where the effective fermion interaction has 
been eliminated by the introduction of auxiliary bosonic flelds. 



^^We note that the renormalization group investigation in [78| of a Unear sigma model in 4 — e 
dimensions misses this property. 
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FIG. 11. The zero temperature effective potential U (in MeV ) as a function of cj = (p/2)^' ^ for 
different chemical potentials. One observes two degenerate minima for a critical chemical potential 
/ic/Afg = 1.025 corresponding to a first order phase transition at which two phases have equal 
pressure and can coexist (Mq = 316.2 MeV). 



The extent to which the transition from nuclear matter to quark matter in QCD differs 
from the transition of a quark gas to quark matter in the quark meson model or NJL- 
type models has to be clarified by further investigations. At the phase transition, the quark 
number density in the symmetric phase n]J^ = 0.52Mq turns out in this model to be not much 
larger than nuclear matter density, n^/^^^ = 152 MeV. Here the open problems are related to 
the description of the low density phase rather than the high density phase. Quite generally, 
the inclusion of nucleon degrees of freedom for the description of the low density phase shifts 
the transition to a larger chemical potential and larger baryon number density for both the 
nuclear and the quark matter phases |]79|. An increase in /i^ results also from the inclusion of 
vector mesons. This would not change the topology of the phase diagram inferred from our 
model. Nevertheless, the topology of the QCD phase diagram is also closely connected with 
the still unsettled question of the order of the high temperature (/i = 0) transition for three 
flavors with realistic quark masses. If the strange quark mass is too small, or if the axial 
f/(l) symmetry is effectively restored in the vicinity of the transition, then one may have a 
first order transition at high T which is driven by fiuctuations [Q (see also the discussion in 
section |VE|) . 

Finally, the low-density first order transition from a gas of nucleons or the vacuum to 
nuclear matter (nuclear gas-liquid transition) can only be understood if the low-momentum 
fermionic degrees of freedom are described by nucleons rather that quarks ||79|]. The low- 
density branch of figure |1^ cannot be carried over to QCD. We emphasize that nucleon 
degrees of freedom can be included in the framework of nonperturbative fiow equations. 
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FIG. 12. The plot shows Uq , where nq denotes the quark number density as a function of /U 
in units of the effective constituent quark mass {Mq = 316.2 MeV). 



E. The tricritical point 



Already from the qualitative information about the order of the chiral symmetry restoring 
transition (a) as a function of temperature T for vanishing chemical potential /i, and (b) as 
a function of /i for T = 0, one can argue that there exists a tricritical point with long-range 
correlations in the phase diagram. The physics around this point is governed by universality 
and may allow for distinctive signatures in heavy ion collisions [^,0,|l^]. We will concentrate 
first on the "robust" universal properties of this point. 

In section [IV B| we find that at /i = 0, the chiral phase transition as a function of 
increasing temperature is second order for two massless quarks. In the chiral limit (ipip) 
vanishes identically in the high temperature phase whereas it is nonzero for T < T^.. A 
singularity separates the two phases which can also be seen from the fact that the function 
{ipip) = cannot be continued to (ipip) ¥" analytically. In the vicinity of this singularity the 
physics is governed by universality. There are only two independent critical exponents near 
this ordinary (bi) critical point. There is no reason to expect a small chemical potential to 
change this result, since this introduces no new massless degrees of freedom in the effective 
three dimensional theory which describes the long wavelength modes near T = Tc. Since 
the order parameter (ipip) is identically zero on the one side of the transition, the line of 
phase transitions emanating from the point T = Tc, fi = cannot end in the T/x-plane. On 
the other hand, in the previous section we find that chiral symmetry restoration at T = 
proceeds via a first order transition within the considered models. Therefore, the minimal 
possibility is that the line of second order transitions coming from T = Tc, fi = turns into 



a first order transition at a specific value T 
where this occurs is a tricritical point. 



Ttc, /i = /itc- The point in the phase diagram 
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Let us first consider the physics around this point in terms of the effective Landau- 
Ginsburg theory for the long wavelength modes ~ {ipip)- In the vicinity of T = Ttc, 
/x = /itc this requires a 0^ potential which has the form 

t/(0, 0; /., T) = f/(0, 0; /z.^, T.J + ^^0^ + ^^0^ + ^^0^ - /.0 , (96) 

where the coefficient h of the linear term is proportional to the current quark mass. The 
coefficients a and h are both zero at the tri critical point, and barring accidental cancellations 
both will be linear in both (T — Ttc) and (/i — /itc)- The [i and T dependence of c is not 
important, as c does not vanish at the tricritical point, and it is convenient to set c = 1. It 
is easy to verify that, near the tricritical point, the line of second order transitions is given 
by a = 0, 6 > and the line of first order transitions is given by a = 36^/16, 6 < 0. Minima 
of VL are described by the scaling form 



h = ^l\-l + -o+A- (97) 




From this, we read off the exponents (5 = 5, 1/ (3 = 4, and 0t//? = 2 where 0t is called 
the crossover exponent because tricritical (as opposed to first or second order) scaling is 
observed for b < a'^'K For more details see ref. [^. The critical fluctuations are described 
by an effectively three dimensional theory (cf. sect. IV A ) and the 0^ coupling in (^) is di- 



mensionless, becoming a marginal operator. This explains why the upper critical dimension 
for the 0^ theory is c? = 3, and the behavior in the vicinity of the tricritical point is correctly 
described by mean field exponents up to logarithmic corrections [^. So if we trust the 
qualitative feature that the transition is second order at high temperatures and first order 
at low temperatures, then QCD with two flavors will have a tricritical point in the same 
universality class as that in our model, with the critical exponents given above. 

Physics away from the chiral limit. What, then, happens in this region of temperature 
and chemical potential in the presence of a small quark mass? We have seen that a nonzero 
quark mass does have a qualitative effect on the second order transition which occurs at 
temperatures above Ttc: The 0(4) transition becomes a smooth crossover. A small quark 
mass cannot eliminate the first order transition below T^. Therefore, whereas we previously 
had a line of first order transitions and a line of second order transitions meeting at a 
tricritical point, with m 7^ we now have a line of first order transitions ending at an 
ordinary critical point. The situation is precisely analogous to critical opalescence in a 
liquid-gas system. At this critical point, one degree of freedom (that associated with the 
magnitude of 0) becomes massless, while the pion degrees of freedom are massive since chiral 
symmetry is explicitly broken. Therefore, this transition is in the same universality class as 
the three-dimensional Ising model. 

With the nonperturbative methods presented in this lecture we can obtain the equation 
of state near the critical endpoint, in the same way as we obtained the universal properties 
of the high temperature (/x = 0) transition in section [IV C|. If we are just interested in 



the universal properties we can directly study a three dimensional model with a single 
component field. The flow equation for the effective potential in the vicinity of the critical 



endpoint corresponds to (|5^) if the fermionic contribution ~ A^^ and the contribution from 



the three pions are neglected. The only contributions comes then from the a-resonance. 
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This flow equation for d = 3 lias been solved numerically in ref. [0 and detailed results and 
comparison with other methods can be found there. 



From our studies of the chiral phase transition at nonzero temperature in section [IV A 



we 

have observed the possibility of a second order transition, with infinite correlation lengths 
in an unphysical world in which there are two massless quarks. It is exciting to realize 
that if the finite density transition is first order at zero temperature, as in the models 
we have considered, then there is a tricritical point in the chiral limit which becomes an 
Ising second order phase transition in a world with chiral symmetry explicitly broken. In 
a sufficiently energetic heavy ion collision, one may create conditions in approximate local 
thermal equilibrium in the phase in which spontaneous chiral symmetry breaking is lost. 
Depending on the initial density and temperature, when this plasma expands and cools it 
will traverse the phase transition at different points in the (/i, T) plane. Our results suggest 
that in heavy ion collisions in which the chiral symmetry breaking transition is traversed at 
baryon densities which are not too high and not too low, a very long correlation length in 
the a channel and critical slowing down may be manifest even though the pion is massive. 
We observe the tricritical point (critical endpoint) to emerge from the qualitative result 
of a second order transition (crossover) in one region of phase space (high T, small /i) and 
a first order transition for the same order parameter in another region (high fi, low T). 
The universal properties of this point can be determined quantitatively. The location of the 
tricritical point in the phase diagram is a nonuniversal property, and the crude knowledge 
of the critical chemical potential for the high density, zero temperature transition makes a 
precise estimate difficult (cf. section |VD| ). The values for T^ and fitc have been estimated 



within a NJL-type model in a mean field approximation to be about /itc — 200 MeV and 
T = Ttc — 100 MeV |0, which agrees rather well with an estimate using a random matrix 
model ig. 



Finally, the question of whether the tricritical point is realized in QCD is closely con- 
nected with the question of the order of the high temperature (/i = 0) transition for realistic 
strange quark masses. If the strange quark mass turns out to be too small then we may 
have a first order transition at high T which is driven by fluctuations |@]|3- As a conse- 
quence, for a particular value of the strange quark mass, rris = m*, one expects the presence 
of a tricritical point even for fi = ||9|- [rTl , |58| . In this case and also for rris < m*, a line 
of first order transitions connects the T and /i axes. Endpoints would then only occur if 
the high temperature, low density region and the low temperature, high density region are 
disconnected. 



^^The same situation arises if the axial U{1) symmetry is effectively restored about the transition 
temperature. 
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F. Color superconductivity 



Chiral symmetry breaking in vacuum is due to a quark-ant iquark pairing with zero net 
momentum. As the density grows, more and more low momentum states are excluded from 
pairing and the chiral condensate is suppressed. We have verified this behavior in section |VD| 
within the effective quark meson model. On the other hand, there are attractive quark-quark 
interactions. In contrast to quark-antiquark pairing there is little cost in free energy for 
correlated pairs of quarks (or antiquarks) near the Fermi surface. This situation is visualized 
in figure |13|. (Figure taken from ||8^.) As in the case of ordinary BCS superconductivity in 




FIG. 13. The figure shows a comparison of quark-antiquark (a) and quark-quark pairing (b). 
Arrows indicate typical pairs, having net zero momentum. 



electron systems, an arbitrarily weak attractive interaction between quarks renders the quark 
Fermi surface unstable and leads to the formation of a condensate. Since pairs of quarks 
cannot be color singlets, a diquark condensate breaks color symmetry. It was originally 
suggested by Bailin and Love (see also [^]) that QCD at very high density behaves as a 
"color superconductor" : Cooper pairs of quarks condense in an attractive channel opening 
up a gap at the Fermi surface. Recent mean field analyses of NJL-type models indicate 
that quark pair condensation does indeed occur and that the gaps are phenomenologically 
significant (of order 100 MeV) at densities of about a few times nuclear matter density [Q . 
Different attractive channels may lead to a possible (simultaneous) formation of condensates. 
In particular, the vacuum of QCD already has a condensate of a quark-ant i- quark pair. An 
important ingredient for the understanding of the high density phase structure is the notion 
of competing condensates. One expects, that the breaking of color symmetry due to a {ipip) 
condensate is suppressed by the presence of a chiral condensate. Likewise, chiral symmetry 
restoration may be induced at lower densities by the presence of a color superconductor 
condensate. This behavior has been demonstrated explicitly within an NJL-type model 
and led to a first model calculation of the phase diagram 0, similar to the one sketched 
in the figure of section |I[ At low temperatures, chiral symmetry restoration occurs via a 
first order transition between a phase with low (or zero) baryon density and a high density 
color superconducting phase P, p4| , ^ . Color superconductivity is found in the high density 
phase for temperatures less than of order tens to 100 MeV. The strong competition between 
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chiral and superconducting condensates also simplifies the discussion: The behavior of the 
chiral condensate within the low density phase is well described without taking into account 
diquark condensation once the phase boundary is known. We implicitly exploited this fact 



in section VD. 



A systematic calculation in perturbation theory is only possible at asymptotically high 
densities, where the strong gauge coupling at the Fermi energy is small. The corresponding 
renomalization of quark operators in the vicinity of the Fermi surface has been studied 
in [Q. The method using nonperturbative flow equations, which were discussed in these 
lectures, may help to give some insight about these questions at lower densities. Since this 
remains to be done I will content myself here with some general aspects and an outlook. 

The explicit form of the favored diquark condensate crucially determines the symmetry 
properties of the ground state at high density. Diquark operators have the structure 



i.Ts 



V^'s 



Dirac ^ Color ^ Flavor 



^ 



(98) 



where the combination of Dirac, color and flavor matrices have to be antisymmetric according 
to the Pauli principle. It is instructive to do some classification. Let us divide matrices 
in symmetric and anti-symmetric ones, where we consider the case of vanishing angular 
momentum 



-'Dirac 



"'Dirac 



color 



yS 
color 



-'flavor 



■'flavor 



C75, C, C7^75 



-^7) "-^5) A2 > = S A 
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(A 



A\ 



jjk 



—te. 



ijk 



Ao = Vol) A3, Ag, Ai, A4, Ag 






r2 



^2. 



«j 



le-ij 



1 
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(99) 



Here "A" and "S" denotes antisymmetric and symmetric, respectively; the r's are the Pauli 
matrices; C is the charge conjugation matrix, E^j^j. are the three antisymmetric and E^^ior 
the six symmetric Gell-Mann matrices. This split corresponds to the fact that the two 
quarks can be comined to form a color antitriplet and a symmetric color 6: Zc®^c = 3c + 6c. 
Consider first the color 3 representation. Overall antisymmetry then requires that combined 
Dirac-flavor has to be symmetric. If we consider the antisymmetric flavor singlet T2 then 
the Dirac matrices must be either C75, corresponding to a Lorentz scalar condensate, or 
C which is a pseudoscalar. The other possible combination involving Cj^'-f^ does not pair 
states near the Fermi surface and is suppressed. Instead it realizes pairing of a state near 
the Fermi surface and one much below with net zero momentum. 

Let us turn to the question if there is indeed an attractive channel between quarks of 
different color. We may content ourselves here by looking at the signs for the scalar channels 
in the effective four quark interaction Ai ([2l|) 
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-M(Pi,P2,P3,P4) = -^A:(-Pi)7Mf i^]i-P3)\\i^tip^hA^] 4iP2)\. (100) 



The curled brackets indicate contractions over spinor indices, i,j, k, I = 1,2,3 are the colour 
indices and a,b = 1,2 the flavour indices of the quarks. Let us neglect for a moment the 
color and flavor structure of the interaction. Using the identity 

+ \{^l'C)^,{^Cl,)8p - \{l'YC)^,{C^'^^,)sp (101) 

we can rewrite the interaction in terms of quark and antiquark bilinears^, respectively 
M{pi,p2,P3,P4) = +{^{-Pi)iC^'^{p4)} {ilJ^{p2)iCtlj{-p3)} 

+ {H-Pih'cij'^ip^)} {ViP2)c^'i^i-P3)} 

+ {H-Pih^YC^p'^iPi)} {V{P2)Cj'j,^p{-P3)} . (102) 

From the signs in (|102|) one may expect repulsion in the pseudoscalar and scalar diquark 
channels (compare also with the attractive scalar and pseudoscalar meson channels in (p^ ). 
However, we neglected color and flavor so far. The appropriate Fierz identity for color 
manifests the separation in antitriplet and color 6 diquark channels 



For flavor we use 



1 ^ 

SabScd = -i^2) ac{r2)db + J2(^!)U^z)db ■ (104) 

^ z=l 



From the different relative sign for the 3c and the 6c channels in ( |10^j| ) we conclude that one 
may expect condensation in the color antitriplet Lorentz scalar and pseudoscalar channels. 
The corresponding color symmetric channels may not be expected to condense. 

We note, however, that the four quark interaction (^) ("dressed one gluon exchange") 



does not break the axial U{1)a of QCD. It has been shown pj] that the instanton induced 
interaction between light quarks properly reflects the chiral symmetry of QCD: axial baryon 
number is broken, while chiral SU{2)l x SU{2)r is respected. (The corresponding f/(l)^- 
breaking term can be written as a determinant (p8|) which has been taken into account 



in the effective quark meson model in section III.) Here we only note that the instanton 



^^In our Euclidean conventions |^] {i^iCiP")'^ = (ip'^iCil^), (ipj^C'ilj'^)'^ = (ip'^Cj^ip). Here f 
denotes the operation of Euclidean reflection in analogy to hermitean conjugation in Minkowskian 
space time. 
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induced interaction leads to an attractive scalar diquark channel but exhibits repulsion in 
the pseudoscalar channel for 3c 0,§ . For two flavors triplet pairing of the form 

mp)C^,i^i{-p)) = K{p^)e'^'ea, (105) 

has been demonstrated to lead to significant gaps of the order 100 MeV for instanton induced 
interactions [^Q- For a nonzero expectation value of ( |105| ) the diquark bilinear has a 
color index i which chooses a direction in color space, say £ = 3, thus breaking color 
symmetry SU{?)) —* SU{2). It leaves all flavor symmetries and, in particular, the chiral 
SU{2)l X SU{2)r intact. One may expect it to be especially favorable because it maintains 
a large symmetry, so that the interaction of a given pair can obtain contributions from other 
pairs with different particle combinations. The above pairing does not involve quarks of the 
third color. It is possible that the remaining ungapped Fermi surface becomes unstable too, 
which is of course color symmetric. An additional condensation of the form 

{iPl{p)Cao,4j't{-p)) = v{p')PiSab (106) 

has been investigated 0. This condensate violates rotation symmetry and the gaps were 
found to be of order several keV at best. 

Breaking of the color SU{3) symmetry generates Goldstone bosons, formally. However, 
since color is a gauge symmetry and the Higgs mechanism operates, the spectrum does not 
contain massless scalars but massive vectors. The number of broken generators is 8 — 3 = 5 
for condensation of the form (|105|) and, therefore, five gluons become massive. We have 
decribed here the color superconducting phase as a Higgs phase. One expects, however, 
that there is a complementary description in which this is a confining phase. Indeed, the 
color superconducting phase for two flavors may be considered as a realization of confinement 
without chiral symmetry breaking. 

The U{1)q of electromagnetism is also spontaneously broken by the condensate. How- 
ever, the photon mixes with the eighth component of the gluon to generate an unbroken 
U{1)q>. Let T3 denotes the diagonal generator of isospin SU{2) and B baryon number then 
Q = T3 + B/2. For a nonvanishing color antitriplet condensate ( |105| ) A ~ e*-''^ one finds 



BA=^A, QA = ^A, A8A = -^A (107) 

and one observes a linear combination of electric charge and color hypercharge Y = X^/y/S 
under which the condensate is neutral: Q' = Q — Y/2, so a modified U{1) persists. (Instead, 
one may also want to refer to baryon number and finds the combination B' = B — Y which 
leaves the condensate invariant.) 

Three quark flavors. The situation is further complicated by the question about the 
number of light quark flavors participating in the condensate. If the strange quark were 
heavy relative to fundamental QCD scales, the idealization of assuming two light flavors 
would be obviously sufficient. The situation becomes more involved since we are interested in 
densities above the transition for which the chiral condensate vanishes. Without spontaneous 
chiral symmetry breaking the relevant strange quark is of the order 100 MeV. Since we are 
interested in chemical potentials larger than 100 MeV it is a quantitative question if an 
approximation with two light flavors is realistic. 
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In the limit of three massless quark flavors a very interesting and compelhng possibihty 
for condensation has been proposed by Alford, Rajagopal and Wilczek 0. Consider a 
generahzation of the two flavor condensate ( |105| ) where the antisymmetric flavor singlet 
~ eab is replaced by the antisymmetric triplet for three flavors ~ eatc Since there are an equal 
number of flavors and colors there is the possibility for a color-flavor structure ~ e^-^^Cab^ = 
61^61 — 6l6i by summing over z. In a slightly more general variant 

{^:{p)Ci,i^i{-p)) = m{p')6l6i + K,{p')6l6i (108) 

has been investigated. The mixed Kronecker matrices are invariant only under matched 
vectorial color/flavor rotations. Ignoring electromagnetism the symmetry of QCD with 
three massless flavors is 5'f7(3)coior x SU{3)l x SU{3)r x Ub{1) (cf. section [11 A| ). The 
SU{3)coioT is a local gauge symmetry while the chiral flavor symmetry is global, as is the 
U{1)b for baryon number. As a consequence of the above matching between color and 
flavor symmetries for a nonvanishing condensate the only remaining symmetry is the global 
diagonal SU{3)co1ot+l+r subgroup. In particular, the gauged color symmetries and the global 
axial flavor symmetries are spontaneously broken. The latter means that chiral symmetry 
is broken by a new mechanism. The breaking of the U{1)b of baryon number leads to 
superfluidity as in liquid helium. All fermions aquire a gap and there are Goldstone modes 
associated with the spontaneous chiral symmetry breaking. The high density phase may 
therefore be in many ways quite similar to low density QCD ||^! 

It is not clear what phase will be realized with physical mass u, d and s quarks. As 
pointed out in |^ , one may think of having the two flavor condensate flrst and color-flavor 



locking at higher density. A flrst investigation may be performed along the lines presented 
in where simultaneous condensation phenomena have been discussed, but flnally more 
sophisticated methods than the mean fleld analysis will be needed to settle these questions. 
An interesting possibility will be the use of truncated nonperturbative flow equations along 
the lines presented in these lectures. 
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